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We perform a systematic study of various versions of massive gravity with and without violation of Lorentz 
symmetry in arbitrary dimension. These theories are well known to possess very unusual properties, un- 
familiar from studies of gauge and Lorentz invariant models. These peculiarities are caused by mixing of 
familiar transverse fields with revived longitudinal and pure gauge (Stueckelberg) fields and are all seen 
already in quadratic approximation. They are all associated with non-trivial dispersion laws, which easily 
allow superluminal propagation, ghosts, tachyons and essential irrationalities. Moreover, coefficients in front 
of emerging modes are small, what makes the theories essentially non-perturbative within a large Vainshtein 
radius. Attempts to get rid of unwanted degrees of freedom by giving them infinite masses lead to DVZ dis- 
continuities in parameter (moduli) space, caused by un-permutability of different limits. Also, the condition 
nigh = oo can not be preserved already in non-trivial gravitational backgrounds and is unstable under any 
other perturbations of linearized gravity. At the same time an a priori healthy model of massive gravity in 
quadratic approximation definitely exists: provided by any mass level of Kaluza-Klein tower. It bypasses the 
problems because gravity field is mixed with other fields, and this explains why such mixing helps in other 
models. At the same time this can imply that the really healthy massive gravity can still require infinite 
number of extra fields beyond quadratic approximation. 
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1 Introduction 



Renewed interest to massive gravity [1, 2, 3, 4, 5, 6] and its further modifications, involving mixing with extra 
light fields and/or tiny violation of Lorentz symmetry is dictated by problems of nowadays cosmology, caused 
by spectacular advances of observation astronomy. These days massive gravity is one of the so-far-desperate 
attempts to cook up a theory which naturally explains the phenomenology of hidden energy, which is currently 
thought to account for over one half of the energy density of the Universe. 

Whatever is their relevance for phenomenological purposes, the problems of massive gravity are of deepest 
theoretical interest. Abandoning Lorentz invariance one actually opens a Pandora box of hidden structures, un- 
observable in the massless case. They include a variety of dispersion (spectral) relations for different components 
of the graviton field, almost as rich as in solid state physics, with non-quadratic dispersion laws, superluminal 
propagation, emergency of non-trivial spatial structures (wave densities) etc. 

This paper (tightly connected with [6]) is an attempt to understand in the most primitive linear- algebra 
terms the puzzling properties of massive gravity [7, 8, 9] and of the whole new world arising after the violation 
of Lorentz and gauge symmetries, which was discovered in [10, 11, 12] and nicely reviewed recently in [13] (see 
also [14]). Since reasons for this strange behavior are not the main concern for all these papers, which are more 
interested in enumeration of diflFerent models and their relevance cosmological applications, our goal is to make 
a step in this direction. Actually we are going to perform analysis in the old-fashioned style of [15] and of [16], 
where the similar puzzles of topological massive photodynamics [17] were addressed and resolved. A posteriori 
it looks quite similar in spirit to the original presentations in [7] , and also includes a direct generalization of the 
original DVZ-approach to the case of Lorentz-non- invariant theories and/or models with extra scalar fields. It 
is of course very close to the original papers [10, 11, 12, 13], just our accents are different. Our interest to the 
problem was initially motivated by studies of massive graviton radiation [18]-[19] in another class of speculative 
models, related to micro rather than macro world: in the TeV scale gravity [20], where masses are presumably of 
Kaluza-Klein origin and various problems of the massive gravity are supposed to be absent (other manifestations 
of the TeV scale gravity are discussed in [21]- [22]). 

In the present paper we analyze only quadratic approximation to Einstein-Hilbert Lagrangian (linearized 
gravity), thus all efi'ects of field interactions, including Vainshtein radius [8] or Boulware-Deser modes [9] and 
superluminal effects [23] in curved backgrounds are beyond the scope of this text. As already explained in [13] 
they are in fact intimately related to peculiar properties manifest in quadratic approximation, in particular to 
the DVZ discontinuity [7] (see also [24]). 

Given a quadratic action, one immediately obtains the Born interaction between currents: 



After Fourier transform the entries of Kjj are quadratic polynomials in the frequency co and space momentum 
k (with some mass terms added), and 

detK = '[[Xaiu>,k) (12) 

a 

In the Lorcntz-invariant case the fc-dependence is of course related to w-dependence and Aa(a;, k) — > Xa{k^) with 
fc^ = -w^ + k"^ so that 

deti^ ^^2_A^(fc) ^•'^^ 

however, Lorentz violation leads to more sophisticated denominators in (1.3). Coming back to the Born inter- 
action, one can rewrite it as 

JK-'J^Y. "t-'y-^' (1.4) 

with some "structure constants" a{k). 

The problem of dispersion relations is basically that of the eigenvalues of K{k): roughly, uj = £(|fc|) is 
a condition that some eigenvalue A(fc) — 0. However, this "obvious" statement requires a more accurate 
formulation. The point is that K is actually a quadratic form, not an operator, what means that it can always 
be brought to the canonical form with only ±1 and at diagonal, thus leaving no room to quantities like A(fc). 
Still, this "equally obvious" counter-statement is also partly misleading, because we are interested not in an 
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isolated quadratic form, but in a family of those, defined over M. This means that the sets of ±1 and can 
change as we move along A4, and degeneracy degree of quadratic form K{k) can change. Of course, this degree 
(a number of O's at diagonal) is an integer and changes abruptly and thus is not a very nice quantity. A 
desire to make it smooth brings us back a concept of X{k). However, in order to introduce A(A;) one needs an 
additional structure, for example, a metric in the space of fields. 

In application to our needs one can introduce "eigenvalues" A(A:) as follows: consider instead of 11 = J^J 
a more general quantity 

n(A|fc) = Jy^jJ (1-5) 

Then as a function of A it can be represented as a sum of contributions of different poles: 

be 
Aa 



n(A|fc) = (1-6) 



then Xa{k) are exactly the "eigenvalues" that we are interested in, and our original 

a,b,c 

The only thing that one should keep in mind is that this decomposition depends on the choice of additional 
matrix (metric) /, which can be chosen in different ways, in particular, its normalization can in principle depend 
on the point of A4 . We shall actually assume that it does not, and clearly the physical properties do not depend 
on this choice, however, concrete expressions for Xa{k) do. It is important, that the dispersion relations - the 
zeroes of Aa (fc) - are independent of I. 

Introduction of I is also important from another point of view. To be well-defined, the Lorentzian partition 
function requires a distinction between the retarded and advanced correlators (Green functions) , which is usually 
introduced by adding an infinitesimal imaginary term to the kinetic matrix K: the celebrated ie in the Feynman 
propagator. However, in the case of kinetic matrix this is not just ie, it is rather ielp with some particular 
matrix Ip. If we identify our I with Ip, then the dispersion relations are actually 

Aa(fc) = ie (1.8) 

what implies that Xa{k) is, in fact, very different from — Aa(fc), and this is related to the important concept of 
ghosts. 

When we have a family of theories, like massive gravities with different masses, we actually have the set 
of quantities Xa{k) and Q;^^(fc) "hanging" over each point of parameter (moduli) space and one is interested in 
the change of this structure when we move around in the moduli space. What happens, different Aa can cross 
or merge, they can also go away to infinity, even more interesting are the properties of a's. At some points 
of the moduli space the symmetry of the underlying theory is enhanced, and we obtain a singularity, where 
limits along different directions do not coincide (this is exactly the reason for DVZ "discontinuity"), so that 
such points should actually be blown up to resolve the singularity. All this is a typical string-theory subject 
[25] , it is amusing that this standard set of questions unavoidably arises in the study of such seemingly innocent 
subject as linearized massive gravity... 

Schematically we consider different intermediate particles (components of the gravity field) , which contribute 
in different channels in (1.4), or, putting this differently, just diagonalize the coefficients a w.r.t. current indices 
so that the sum (1.4) becomes a sum over channels and particles. In general we are interested in a variety of 
channels and in contribution of different species to each of this channel. The whole pattern is characterized by 
the following data (see also [6]): 

Species: dispersion laws. Sometimes the dispersion law is simple. 



\lc^P + M^ (1-9) 



however, the coefficients and are of importance. The dispersion law with > 1 describe superluminals 
and with < tachyons. The superluminals always travel faster than light and can violate naive causality, 
[26], [23]. They are sometimes also called tachyons in literature. But physically, they are very different from 
tachyons, which are signals of instabilities and do not violate causality (allow simultaneous but uncorrelated 
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and causally independent development of instabilities at space-like intervals). In fact we shall see that there are 
more sophisticated dispersion laws 

a; = A(|fc|) ^±V^c2fc2 + M2 (1-10) 

however their complexity does not increase with increase of the dimension d: relations between u) and are 
at most quartic: 

a;2 = sjp^^)±^fp^) (1-11) 

In the case of this more complicated dispersion law, we define as position of the pole in uP' at zero spatial 
momentum, and define the tachyon mass square as a real-valued solution to the equation K{k'^) = 0. 

Residues: coefficients controlling the contribution of the specie a to the channel h. It is important 
to distinguish if ghost contributions appear in physical (say, space-time transverse for conserved currents) or 
unphysical channels (corresponding to sources of pure-gauge species). 

Discontinuities appear when some — > oo. Accurate formulation of the problem is that we look at the 
interaction in a given channel at large distances, but not as large as min[M^^), so that all contributing species 
still look like massless. However, if some Mqq = oo there is simply no such region and we observe a jump 
between "long-distance" interactions for Ma^ = oo and Ma^ = 0. In other words, switching on a tiny mass scale 
is not obligatory a small eff'ect if for some specie this tiny scale is multiplied by infinity, constructed from other 
parameters (like {A/B — in the Pauli-Fierz case below). 

In the fully-comfortable theory there are no tachyons, ghosts and super luminals. This, however, is rarely 
achievable in theories of massive spins, greater than 1, if one decides to abandon gauge invariance and give 
masses by explicit rather than spontaneous breaking of gauge invariance. 

In the paper, we analyze the problem in a somewhat unusual way. Instead of defining the normal modes 
by passing to Hamiltonian formalism we directly diagonalize the kinetic matrices Kjj and Kf in momentum 
representation. The normal modes defined in the first way are sometimes very convenient to deal with, especially 
in the case of massless theories, since despite they are introduced in explicitly Lorentz non-invariant way, the 
longitudinal modes are very distinguished in this case, they carry a lot of physical information, in particular 
easily distinguish between propagating and non-propagating modes (without throwing the latter away from 
the spectrum as one often does by imposing constraints in the Hamiltonian formalism what in fact makes 
deformations to adjacent points in the moduli space problematic, see [16] for initial criticism of the standard 
approach). On contrary, in the massive theories it is more convenient to use the Lorentz invariant normal modes 
(the both types of modes coincide in the rest frame). 

In what follows we first start from a simplest example of electrodynamics in order to illustrate the mode 
based approach and then continue in s.4 with massive and massless gravity. One of our purposes in section 4 
is to justify the above-mentioned mechanism of DVZ discontinuity. Then in s.5 we consider generalizations to 
Lorentz- violated gravity, where the main novelty is occurrence of quasiparticles with no n- trivial dispersion laws, 
not very familiar in elementary particle physics. Finally in s.6 we turn to massive gravity mixed with some extra 
particles. Kaluza-Klein massive gravitons belong to this class, what a priori explains why addition of mixings 
can produce healthful theories of massive gravity. For other fashionable models of this kind see [11, 4]. The last 
section contains a discussion of various physical consequences of the behaviour obtained in previous sections. 
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2 A warm-up example: massive photodynamics 



In this section we consider the case of photodynamics which is much simpler than the gravity and, hence, we 
use it to ilhistratc the approach of the paper. We look at various patterns of adding masses, including those 
breaking the Lorentz invariance. 



2.1 Generalities 

Photodynamics is the theory with quadratic action 

J (^-^F^.^F'^'' -M^Al + J^A''y''x = j (A^{-k)K^„{k)A''{k)+A^'{-k)Ji,{k))d''k (2.1) 
Our immediate task is to: 

- enumerate different modes, contained in the field A^j^, which propagate through the space- time indepen- 
dently, without mixing, and identify their properties, 

- find the sources of these modes, 

- decompose the Born interaction between the sources into contributions of different modes. 

In this way we can express various properties of interaction, mediated by our theory, through the properties 
of individual modes and in this way identify the origins of particular types of unusual behavior. 

Of course, formal realization of this "program" is nothing but an elementary linear algebra exercise with the 
kinetic matrix K^^{k) in momentum representation, which in the case of photodynamics is simply an ordinary 
symmetric dx d matrix: 

K^oc = k^kc, - {k^ + M^)r)^^ (2.2) 



It can be easily diagonalized: 



where the d eigenvectors are: 



a=l 



gauging scalar v'^ = qk^ 

longitudinal vector v\^^ (2.4) 

transverse vector v^, i = 1, . . . ,d — 2 

"Transverse" and "longitudinal" refer to space rather than space-time vectors. All these components are well 
defined at 7^ and the splitting exhibits a smooth limit in the massless limit 0. 

The gauge degree of freedom, A/j^ = qk^ is scalar, it does not mix with the other d — 1 degrees of freedom, 

k^'K^^ = -M^kc, (2.5) 

(it is a particular eigenvector of K), but it has a non-trivial Lagrangian and even a kinetic term for the 
Stueckelberg field q{x), whenever ^ and gauge invariance is broken. The d — 1 "physical" degrees of 
freedom in their turn split into l+{d — 2) components - longitudinal and transverse photons with different 
eigenvalues and different properties. 

Note that the very definition of normal modes is not Lorentz invariant (even if the theory is): they solve 
an equation K^^v" = A?;'' and not K^aV°' = This is important for making Lagrangian diagonal, when 
expressed through the normal mode - this follows from orthogonality of matrix eigenvectors (not that in the 
case of Euclidean signature there would be no difference - but the difference between propagating and non- 
propagating modes is not seen there). 
Born interaction, 

' J''{-k)P^^{k)r{k)d'^k (2.6) 



/• 



is defined in terms of the propagator - inverse of kinetic matrix - which is well defined for ^ 

p _ _ + -M^ (2.7) 
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Our notation will be as follows: 

Minkowski metric 77,^,^ = diag(— 1, +1, +1, . . . , +1), 

Space-time momentum = (w, fc|| , 0, . . . , 0), = (— w, 0, . . . , 0), spatial momentum will also be often 

denoted by k and frequency w - by fco- 

The Lorcntz square of space-time momentum is fc^ = — o;^ + fc^ = —to'^ + fc^. 

Space time indices are denoted by Greek characters: = 0,1, . . . , d— 1, space indices - by Latin characters 
i,j = 1, . . . ,d — 1, transverse spatial indices (in coordinate where k is directed along the first axis) - by 
a,b = 2, . . . ,d — 1. 



2.2 Massless photon 

This is the ordinary massless photodynamics. 
Kinetic matrix is: 



/ + A;^ u>k 



II 



k^ 









y 

Eigenfunctions and eigenvalues: 



II 



V / 



Ag = 



/ V 



/C|| 










ujk\\ 
















2 7 2 

LO — k^^ . 

















gauging photon = Stueckelberg scalar 



longitudinal photon 



(2.8) 



(2.9) 



V / 
/ \ 

1 



V / 



/ \ 




V 1 J 



Xj_ =iO^ - kn 



transverse photon 



d—2 polarizations 

Thus d—2 transverse photons are just ordinary massless particles with the normal kinetic term, described by 
the eigenvalue — fc^ , which vanishes on-shell, when a; = fc| | . This dispersion law allows non- vanishing values 
of u), what implies that once emitted such particles can propagate by "themselves". 

Pure-gauge photon (gauging scalar or Stueckelberg scalar) has vanishing eigenvalue, this means that it 
completely drops out of the action: it " does not exist" , or is ultralocal: just equals to its source. 
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Longitudinal photon is a non-trivial field, but the eorresponding eigenvalue is positively defined and vanishes 
(i.e. is on shell) at a single point lu ^ =0: dispersion law is simply uj = 0. This means that this field can 

not exist " by itself , it is fully driven by the source. Still, eigenvalue depends of lu and k, what means that even 
for a point-like source the field can be spread in time and space. 
ExpEinding the gauge field in different sorts of photons. 



and substituting it into the action / [Ffj^^F^'^ + A'* J^), we obtain: 

Al{uj'^ + klf + A\{uj'^ - kl) + Ag{-ujjQ + fc|| J||) + A||(fc|| Jo + wJii) + A_L Jj_ 



(2.10) 



(2.11) 



Note that we work with un-normalized eigenvectors: this simplifies some formulas, though makes some other 
like this one for the Lagrangian look a little unusual. 

The coupling AgJ vanishes if is conserved, in this case also the A\\J coupling can be rewritten as 



A||(fc|| Jo + a;J||) 



J|| = II 



w " " k\\ 

so that the action, expressed in terms of the separated variables, becomes 

,2 I ;,.2 



^11 Jo 



(2.12) 



Aliu' + klf+ Aliuj' -kl) + ^i^i- J|| +A^J^ = 



,2 I ^,2^ 



^^^2 _|_ 

= Aliuj^ + klf + Al{w^ -kl)+ ,^ " A||Jo + A^J^ 



(2.13) 



Born interaction can be immediately read from this formula: 



J2 

JnP°"^Ju = CO k\\ = --Tr 



/2 



Of course it can be alternatively obtained by first writing down the propagator P""^ which satisfies 

IC POLV _ CI/ '>'"'>' 



and is equal to 



poiv _ 



fc2 

+ ck°'k'' 



so that the interaction of two conserved currents, such that A;'' J^j^ = 0, is 



;2-fc2 



II 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



In order to obtain (2.14) from this simple expression one should substitute an explicit resolution of conservation 
constraint for the current: 



/ Jo \ 



J 



/ ak\\ \ 



\ J± J 



(2.18) 



where a = # = 



Note that because of the gauge invariance (i.e. the vanishing of eigenvalue Xg) the propagator is not just 

inverse of kinetic matrix, one should exclude the zero mode by putting transverse matrix at the r.h.s. of (2.15). 
Instead propagator (2.16) contains unspecified coefficient c, which drops out from coupling to conserved current. 
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2.3 3d topologically mcissive photodynamics [17, 16] 

Consideration of the massive photodynamics we start with the celebrated "intermediate" example on the way 

from massless to massive photodynamics, when photon gets a mass, but gauge invariance is still unbroken and 
Stueckelberg fields do not show up in the action. Despite propagating photon is massive there is a pole at 
vanishing momentum in the interaction of currents. Remarkably, this does not contradict unitarity because of 
existence of a single propagating massless mode (not a particle) [16], Since it is not a particle, the long-range 
interaction is topological, namely in this case it is just an Aharonov-Bohm interaction. 
Lagrangian: 



(2.19) 



Kinetic 3x3 matrix: 



LOkll 



\ -iMk\\ 



II 



/e|| 



-iMlo 



iMk\\ \ 
iMto 
-e J 



/e|| 



\ -iMk\\ 



UJ 



iMoj 

-iMto - ku j 



(2.20) 



The matrix is not fully symmetric, since the required symmetry property is K^i,(}i) = Kni^^—k). Note also the 
appearance of i = and that 

Eigenfunctions and eigenvalues: 



ej^ because of the Minkowski metric. 



( -u\ 



\ / 



Ag = Stueckelberg scalar 



«ll 

UJ 



/r2 



/ ku \ 



oj^ +r 



■ U! 



photon 



(2.21) 



photon 



where = fc^ + M^{u)^ + kf^). 

Note that the eigenvalues cj^ ± r are quite complicated, however, the dispersion relation ui^ 
equivalent to the standard one: 



2 7 2 
UJ = fcji 



i.e. non-vanishing eigenvalues can be rewritten as 

2^ 

UJ ±r = 



M2)(u;2 



so that no irrationalities show up in denominators in the propagator 



p = 



+ 



r is 



(2.22) 



(2.23) 



(2.24) 



For conserved currents with fc^ J** = 



2 2iMJ_LJ|| 
J I 



UJ 



+-Jl1l = 



U) 



1 



- fc2 _ M2 



2 UMJ^Jq 
J I ; 



t2 
1,2 -^0 



(2.25) 



9 



Note that inseparable products like fc^ (fc^ + Af^) appear in denominators: this will be a typical feature of 
all massive gauge theories, which does not depend - as we already see, on whether gauge invariance is preserved 
or not. 

Not only the modes are inseparable, there is a pole at p|| — 0, despite particles are massive. This pole is 
indeed a long-range interaction and it is a question how such interaction can occur in the theory of massive 
particles. The answer is [16] that the spectrum of the theory is not exhausted by massive particles, there is an 
additional single propagating mode without a gap. This is a single mode, not a particle with the dispersion law 
w = consistent with the Lorentz invariance (see more examples of such dispersion laws in the gravity case 
with broken Lorentz invariance, s.5, especially, s.5.6). Because of this the long-interaction which it describes 
can not transfer space momentum and is especially simple: it is a topological Aharonov-Bohm interaction, but 
it exists, its effects are observable and it can cause infrared divergencies in scattering cross-sections, just as the 
ordinary long-range Coulomb interaction does in low dimensions. 

As we saw, the irrationality, r is gone from denominators in the current-current interactions, moreover, it is 
not seen in the numerators. This, however, is an illusion, irrationalities are excluded at the price of considering 
non-diagonal interactions. If the interaction is diagonalized, i.e. written in terms of independent modes 
(polarizations) and v^, the irrational r is explicitly present in the formulas. This is again a standard 
feature of massive gauge theories. 

In fact, and describe a left-polarized massive photon with two massive degrees of freedom. There is 
an extra degree of freedom (not a particle), responsible for the long-range Coulomb interaction [16]. 



2.4 Massive photon 

This is the theory with the action (2.1)." Here (and only here) we are going to consider for illustrative purposes 
the kinetic matrix not as in (2.2) with both lower indices, but with one lower and one upper. This leads to 
Lorentzian eigenvalues, in contrast to the Euclidean ones obtained from (2.2). This issue is discussed in the next 
section, one can also find the comparison of the Euclidean and Lorentzian eigenvalues in the massive photon 
case in [6]. 

Thus, now the Kinetic matrix is 



u)kii 



—u)k\\ 




-(fc2 + M2) / 



ujk\\ 



-ujk\\ 



\ 



uj^-(kl+M^)) 



(2.26) 



"Since we wish to avoid higher derivatives for Stueckelberg fields (not because they are bad, higher derivatives theory suffers 
from ghosts at most, and even this is not unavoidable [31], sec [36] for a recent summary and a list of references, but simply to 
somehow restrict our moduli space), we do not consider another popular model of massive photon, with idfiAt^)^ term (ironically, 
this was exactly the model analyzed by Stueckelberg in [37]). 
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Its eigenvectors and eigenvalues are 



/ 



ku 



\ 



\ j 



the former Stueckelberg scalar 







^ _ ^2 _ _ the former longitudinal photon 

II II (a scalar) 



(2.27) 



/ \ 



V 1 J 



A_L = — fcn — transverse (d — 2) — vector 



Normalization of modes is such that they polynomial in u and k. 
Gauge field is expanded in different sorts of photons 

Lagrangian, when expressed through the normal modes, becomes diagonal: 
-M^Al + (M2 - a;2 + kfAAf, + {w^ - kf, - M^)Al + J^A^ + " '^""^ 



A;2 -a;2 



-Jok\\ + J\\uj 



(2.28) 



(2.29) 



Accordingly the Born interaction of currents is 

1 M\4 - - Jj) - (Jpg; + J||fc||)2 
4 M2(M2 + fc2 -u;2) 



(2.30) 



If, despite gauge invariance is broken, one keeps the currents conserved, fc^ J'^ = 0, then the Born interaction 
converts into 



1/ J2(/C||-W2) 



. fc2 (M2 + fc2 _ (^2) (7\^2 + j[,2 _ ^2 

Propagator, the inverse of kinetic matrix, is 
/ M2-a;2 



1/ >^||(*ll -'^^) 



w2(M2 + fc2_^2) (M2 + fc2_^2) 



) (2.31) 



u)k\\ 



(w2 -kf,- M'^)M'^ (a;2 - )fc2 _ M2)M2 



ujkw 



I 

fc2 +M2 



(a;2 -kl- M'^)M'^ (^2 - /cj] - M2)M2 



a;2-(A;2+M2) / 



(2.32) 
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and this provides for the Born interaction of currents 
Jo2(/e||-a;2) J2 



4VA;2(M2 + fc2 _^2) 



(M2 + fc2-a;2)) 4(a;2(M2 + fc2 _a;2) 



(M2 + fcj 



-a;2)) 



(2.33) 



what coincides with (2.30) and (2.31). The poles at w = and fc = are spurious, as we already discussed in 
s.2.2. Moreover, in this case there is no pole &i u — k — Q, and thus no long-range interactions as well, as can 

be seen from explicitly Lorentz-invariant formula at the l.h.s. of (2.33). 

Define the propagator as a solution to K^aP°"^ = ^'il ~ ~ Q^)%t-) ^'^ = ~<^^ + ^j^- Then, the smooth 
matching with the massless case (2.15) at a = is provided by the expression for the propagator 



( k'^M'^ - 2M2fc2 + 2fcg 

{-k^ + k'^)M-^{-kl + fc2 + M2) 

kok{M'^ + 2kl) 



(M2 - 2kl)kok 



{-kl + k'^)M-^{-kl + fc2 + M2) 
A;g(2fc2 + M2) 

{-kl + k'^)M'^{-kl + fc2 + M2) (-fc2 + fc2)M2(-fc2 + fc2 + M2) 







-fc2 + fc2+M2 / 



(2.34) 



2.5 Breaking Lorentz invariance 



We assume that the Lorentz invariance SO{d— 1, 1) is broken down to the space-rotation invariance SO(d — 1) 
only by mass terms in the Lagrangian. This means that instead of a single mass term M^A^ = Af2(— + Af) 



there can be two, 



-m2A2 



mjAf, and toq docs not necessarily coincide with mi. If gauge breaking terms 



were added to the kinetic terms also, this would immediately provide fc** terms for the Stucckclbcrg fields in 
the Lagrangian. Without gauge breaking one could only write electric and magnetic terms Fg^ and F^j with 
different coefficients, but this is equivalent to time rescaling and does not really break the Lorentz invariance. 
So, the kinetic matrix is 



/ a;2 + (A;2 + M2) 



V 



u)kii 







ujk\\ 

kf^-ik^+M^) 









\ / k^i+M^ wfcii \ 

wkii - Ml 

V a;2 - (fc^ + Ml) J 



Its eigenvectors and eigenvalues are 

/ 1 

k?, - w'^ -\- + Ml + r 



(2.35) 



2ujk\\ 



A_ = -{uj^ + kl+ - Ml - r) 



fc^ - a;2 + Ml + M| - r 
2wk\\ 





A+ = -{uj^ + kl+ Ml - Ml + r) 



(2.36) 





V 1 / 



uj^-kl- Ml 
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with 

r2 = ^ ^2)2 _ 2(^2 ^ ^2)(^2 _ ^2) ^ (^2 ^ ^2)2 (3.37) 

The gauge field is expanded in different sort of photons: 

A^' := A-v^ + A+v^ + A±v^ (2.38) 
The Lagrangian, expressed through the normal modes is diagonal: 

(P - Qr)A^_ + {P + Qr)Al + {u^ - A:^ - M^)Al + J^A^+ (2.39) 

+ (2Joa;fc|| + J||a;2 - J^k^ - J||(M^ + Mf ) - J^^rj A_ + (^2Joa;fc|| + J^oj'^ - J^kf^ - J||(Mo^ + Mf ) + J^rj A+ 
with 

P = {u^-Miy, Q = Mlk'^^-{2Mf + M^)u^+u^ + oj^k'^^+M^Ml + Mf (2.40) 
and = {uj^+kf,f-2{M^+Mf){ij^-kf,) + {M^+M^f,sothat P^-Q^r^ = 4ij^kf,{oj^M^-k?,Mf-M§M^y. 



The Born interaction of currents is 



Ml 



+ 



^2a;A;|| Jo + (w^ 



II 



Ml - Ml - r 



P-Qr 



+ 



(2ujk\\Jo + (a;2 - fcfj - - Ml + r) J,,)^ 



p + gr 



If the currents are conserved, this converts into 



JliM^uj^-Mlkl) 



JfAM^uj'-Mlk 



a;2 - fc2 - Ml (M2a;2 _ Mffc^ - M^MDk^ _ ^2 _ j^^2 (^^^2^2 _ ^2^2 _ m^MDuj^ 

The propagator, inverse of the kinetic matrix, is 



(2.41) 
(2.42) 



/ 



CO 



Ml 



Lukw 



{M^uj^ - Mlkl - M2) 

LOkw 



{M^uu^ - Mlkl - M2) 



{Mgu^ - Mlkf. - M2) (Mo2a;2 _ Mlkf. - M2) 



V 







(2.43) 



and this provides for the Born interaction of currents: 

J2 j2(M2u;2-M2fc2) 



72 



a;2-(fc2+M2) / 

J|^(Mo2w2_^2^2) 



w2 - fc2 - Ml (M2u;2 - M2A:2 _ M^MDkf^ a;2 _ ^2 _ ^2 (^2^2 _ ^2^2 _ m^MDuj^ 
what coincides with (2.41) and (2.42). 



(2.44) 
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3 Comment on the definition of normal modes 



It is a good moment now to illustrate our view on the normal modes [6]. We see that the physically relevant 
interaction (2.14) contains two very different kinds of structures. In transverse channel interaction is clearly 
mediated by a massless photon: there is a pole whenever this photon is on-shell, = fc^. At the same time 
in the temporal- longitudinal channel the mediator is something very different: the pole of the corresponding 
propagator is at w = fc = and nowhere else, i.e. has (real) codimension two. The seeming codimension-one 
pole at a; = is spurious: as A: — > at w ^ the Jg in the numerator simultaneously tends to zero, so that 
the would-be pole is fully eliminated. Seemingly there is no pole at a; = but fc 7^ 0. These properties of the 
interaction are perfectly encoded in the eigenvalues A: Aj_ = exactly on the mass shell of transverse photon, 
while A|| = only at cj = fc = 0. This is why we consider such definition of A's physically relevant. 

However such definition may look somewhat non-conventional and in fact inconvenient for other purposes. In 
particular it is based on an explicitly Lorentz-non-invariant definition of eigenvectors. An alternative Lorentz- 
invariant definition, however, provides anothcir value of A|| ^ uP' — fc^ and can be made consistent with the 
result for Born interaction only within a sophisticated concept of a "propagating, but decoupling'' longitudinal 
photon. Our approach rather treats longitudinal photon as totally "non-propagating" - in perfect accordance 
with (2.14). 

Technically the difference is as follows. We consider if^^ as an ordinary symmetric matrix and formally 
diagonalize it by orthogonal transformations, i.e. by raising indices with the help of Euclidean , instead of 
rf". In other words we diagonalize the symmetric matrix 



instead of the asymmetric one 



,2 



The pure gauge (gauging) mode with the vanishing eigenvalue is of course the same in both cases, = = 

() , however for the longitudinal mode one gets ^^fi = ( " ) with A|| = lj'^ + k?, instead of the usual 
k J \ J 



= I " j witn /\|| = — /C|| 



with Am = w'^ — kf:. Therefore, our non-covariant modes are orthogonal in the ordinary linear- 



algebra sense, i.e. w.r.t. Euclidean metric (5^,^ instead of Minkowskian ry^^, while the Lorentzian eigenmodes 
are orthogonal w.r.t. 7]^^, i-e. w.r.t. the group SO{d— 1, 1) instead of SO{d). Accordingly, the quadratic-form 
kinetic matrix is 



HI 

in terms of our Euclidean normal mode, and our two eigenmodes, used in (2.10), are just S'^^v]} and e'^^v]}. 
Because of the different A|| our normal mode is clearly non-propagating, while conventional longitudinal photon 
does propagate, just "decouples". As we shall see in the following sections, when mass is introduced, gauge 
invariance broken and the gauging (Stueckelberg) mode revived, the difference gets even more pronounced: 
our normal modes are the mixtures of Stueckelberg and longitudinal modes, one propagating another non- 
propagating, but both with non-trivial eigenvalues (2.27), while in the Lorentz-covariant definitions Stueckelberg 
mode remains non-propagating with the eigenvalue , but longitudinal photon remains propagating, acquires 
the standard dispersion relation w-^ = fc^ + - the same as transverse modes - just it "no longer decouples". 
Thus the standard view is not very helpful in visualizing the actual properties of massive (gauge violating) 
theory with its sophisticated propagators and Born interactions. 

There is a clear physical reason to deal with the Euclidean eigenvalues (see also [6]). Indeed, let us consider 
the massive 4-vector field Lagrangian 

-A^(fc2 + rn')A^ = A^K^^A'' (3.45) 

It is ill-defined as contains ghosts, since the time component of the field has the wrong sign of the time-derivative 
term. This is immediately reflected in the corresponding negative derivative d\{uj'^) / duj"^ of one of the Euclidean 
eigenvalues of K^^, while the Lorentz eigenvalues are all positive in this case. 
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To put it differently, one takes care of the sign in the quadratic action when performing the Gaussian 
integration in the path integral. In particular, integrating over A^, one's concern is only on the coefficients in 
front of Aj^, which are exactly the diagonal elements of the Euclidean kinetic operator. 

It deserves emphasizing that this subtle difference between Euclidean and Lorentz eigenvalues is inessential 
when dealing with only the scalar sector (which is of our main interest throughout the paper). Moreover, 
the majority of results below are independent of this difference between Lorentz-invariant and non-invariant 
definition of normal modes. In particular, independent are det K, characteristic equation and its decompositions. 
Hence, dispersion laws, tachyons etc are independent of the choice of modes, and only the ghost content could 
depend. More comments on this issue can be found in [6], especially in Appendices I-II. 

Note that Lorentz covariant modes are easily restored from our formulas below: in the rest frame (i.e. at 
fc = 0) the two choices coincide and Lorentz transformation can be used to obtain the Lorentz covariant modes 
in all other frames. This is not so easy for non-Lorentz-covariant modes and we have to write them down 
explicitly in all frames. However, even this advantage of the Lorentz modes disappears when one works with 
non-Lorentz-invariant theories, as we do in the second half of the paper. 

In fact, in order to give the reader a flavour of difference between the Euclidean and Lorentzian eigenmodes, 
for illustrative purposes, we present in a couple of places calculations for the Lorentzian modes (as we did in 
s.2.4 above). 



15 



4 Lorentz invariant massive gravity in quadratic approximation 



4.1 Generalities 

Einstein-Hilbcrt action J R^d'^x is highly non-hnear in the metric field and describes a pretty sophisticated 
interacting theory. However, all the peculiarities of massive gravity as compared to ordinary general relativity 
show up already at the level of quadratic action for the small deviation h^^j olg^^i, from the Minkowski background 
rjfj^i,. In this approximation it is also easy to introduce a mass perturbation, which breaks gauge (general 
coordinate) invariance and makes graviton field massive. If Lorentz symmetry SO{d -~ 1, 1) is preserved, the 
possible deviations from General Relativity in quadratic approximation are parameterized by two constants: A 
and B. 

Quadratic kinetic term /i'"'/C^i/,a/3/i"^ is: 



(4.1) 



(4.2) 



-2{k^k^ria0 + kakffTj^J^ - {k'^ + A) (^r/^a??i//3 + Vi^aVfiff^ + 2(fc^ + B)r]^^rjai3 

The propagator - inverse of kinetic matrix - can be parameterized as following 

■p/ii/,a/3 = ai {k^ikailisu + k^kprjav + ki^kaVfiiJ. + ki,kfjr]aij^ + 

+a2k^iku'na/3 + askakisTj^u + ail^rj^ariiyff + fli^a'ntJ.p^ + a^Vni^'naff + a^k^k^kak/s 

and Born interaction of two stress-tensors is given by 

^M.,a/3T'^"r"'5 ^ j,^2 ^ ^ {kTk)T + 2aiT^, + asT^ + ae {kTkf 2aiT^, + a^T^ 

(4.3) 

Note that the Born interaction of conserved energy-momentum tensors that involves only 04 and 05 coeffi- 
cients does not exhaust all possible ways to observe gravity interactions - one can just radiate a graviton by 
the energy-momentum tensor, the radiated field being 

h^,. = V^u^c.pT'^^ = 2ai [k^{kT\ + k,{kT)^] + aafc^fc.T + a^rj^^ (kTk) + 204^+ 

(4.4) 

-Va^^^ivT + aek^k,, {kTk) " — > a2k^kvT + 2aiT^„ + a^ri^i,T 

Of course, violation of gauge invariance (general covariance) in transition to massive gravity liberates stress- 
tensor from obligation to be conserved. However, this implies that violations are also present in the matter 
sector. If, as usual in the present-days discussions, we restrict all violations of conventional physics to the 
pure gravity sector, then the properties of matter are not changed and stress tensor T^^" remains conserved. 
Technically, if one wants to break this property, a new current, J'^ — k^T^^^ , will appear in the formulas, which 
should be further split into conserved Jv = Jv — k^Q with Q = -^k^k^T^^" and "improved", though non-local, 
stress tensor 

f^"" = T^'' -^{k^'J'' + k''J^ + k^k''Q) (4.5) 

will be conserved. 

Therefore, (also in order to allow a smooth transition to the masslcss case) one would better consider not an 
inverse of fC, but solution of the equation ICV = E with the space-time transverse r.h.s., explicitly accounting for 
the conservation of the energy-momentum tensor, which is the source of the gravity field. To make discussion 
complete, we introduce an additional parameter a which interpolates between unity (at a = 1) and space-time 
transverse (at a = 0) r.h.s. E: 

V^^^palC^p = a (77^a??../3 + + {l-a) {vlaVlp + vlfjVla) (4-6) 

where the transverse Kronecker symbol is defined as 

Then coefficients in the propagator are: 
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ai = 



02 



as = 



04 



, ak'^ + aA- A 
'2 Afc2(A:2 + A) 



+ A) - rfAS + (rf - 2)(B - yl)fc2) 

{B + aB- A) + yl.B(a - 1) 
fc2(A;2 + ^) _ dAB + {d-2){B- A)k^^'^ 

1 1 



(4.8) 



05 



2fc2 + A 
{A - B)k'^ + AB 



+ A) - dAB + {d-2){B- 

- a){dB -A)+ a{d - 2){2B - A)k'^ + (1 - a)A [{d - 3)A -{d- 4)B] P 
Ak^{k^ + A) _ dAB + {d-2){B- A)k'^^ 

4.2 Normal modes of the gravity field 

After diagonalization of (4.1) we get the following decomposition of the gravity field /i^i/: 

^ = fc^M±I)+ (.-1) + ^ + 4^ = 

srcoiidary space— time trace ('^•^) 



. . space— I iinv I raii,'^\"cr 
massive spin 2 



Stueckelberg vector 



= < 



d{d - 3) 



+ {d-2) + ^ d_-2 + ^ + ^ ^+1 

longitudinal tensor Spatial trace I I spatial-transverse longitudinal secondary 

, spatial-transverse tensor =t?asverse vector ) (. stueckelberg vector Stueckelberg scalar Stueckelberg scalar 

The first line here describes decomposition into irreducible representations of SO{d— 1) in the rest frame, while 

the second line that w.r.t. helicity group SO{d — 2) acting in the space orthogonal (space transverse) to the 
space momentum k. In fact these decompositions remains relevant even if Lorentz invariance is broken down 
to spatial-rotation symmetry SO{d — 1). In gauge invariant theory Stueckelberg fields do not show up in the 
Lagrangian (this also requires the sources to be transverse). When gauge invariance is broken down by mass 
terms, Stueckelberg fields acquire no n- trivial kinetic terms and they can mix with transverse degrees of freedom. 

In arbitrary frame the two transverse vectors: one from the massive graviton another Stueckelberg get mixed. 
Also mixed arc the four scalars. This means that characteristic equation, which defines eigenvalues of kinetic 
matrix, and is an equation of degree d{d + 1) in fco = w and k is actually factorized: 



Char(A) = (A - A,,)^P2(A)^-2g4(A) = (A - A,.)^(A - A+ ^'^-^(A - K.^-^ A^J = 

a=l 



(4.10) 

where and Q4 are polynomials of degree 2 and 4 respectively and all their coefficients as well as \gr are 

quadratic functions of lo and k. In other words, 
Xgr is some bilinear combination of w and fc, 

^tec ^ P2 ± y/pl, where p2 and are respectively quadratic and quartic in uj and fc, 
Ag^^'^'^ are the roots of degree-four polynomial. 

In the rest frame the roots should be grouped in a different way: 

Char(A) = (A - \r)''^^'^{\ - Xvecf-^X - A+ )(A - A" ) =0 (4.11) 

fe=0 
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i.e. at A; = 



Xtecik = 0) = Xgr{k = 0), 

A^f(fc = 0) = Ag.(fc = 0), 
A|f (fc = 0) = A-,(fc = 0) 



(4.12) 



where "spt" and "sSt" label spatial trace ha and spatial Stueckelberg scalar hij = kikjS respectively. The 

remaining two scalars, the space-time trace (stt) h'^ and secondary Stueckelberg scalar h/j^v = k^kua have 
eigenvalues, which are roots of quadratic equation: 



•^^c = 92 ± \/94|fe=o 



(4.13) 



(4.14) 



In gauge invariant theory all the d Stueckelberg fields have vanishing eigenvalues and we get 
Char(A) = A'^(A - A,.)^(A - A+ J<^-2(A - A^f )(A - A^«)|^^ = 
i.e. in this case the two trace eigenvalues are roots of quadratic equation, 

\t^ = t2±Vu\Gi (4-15) 

Actually gauge invariant will be only the massless gravity (where, by the way, transition to the rest frame is 
not a justified operation). The interrelation between modes in various cases can be described by the following 
table: 



rest frame 



normal modes 



gauge invariant (massless) case 



massive graviton 







graviton 



graviton 



vector 



vector 



Stueckelberg [d — l)-vector 



secondary Stueckelberg scalar 







space-time trace 



Stueckelberg vector 
Stueckelberg scalar 
secondary Stueckelberg scalar 
spatial trace 
space-time trace 



Stueckelberg d-vector 



spatial trace 
space-time trace 



It is also easy to describe the eigenvectors. In coordinate system where k = k\\s directed along the first axis 
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the corresponding modes look as follows (the matrix in the upper left corner is formed by directions and 1): 



V 



u 


u 


u 


u 








... 



















hab 












/ 


A 

u 


U 


... V 


\ 

... 








... w 




V 


w 












/ 


a 






\ 




1 










6 










^1 



transverse traceless graviton 



two transverse (d — 2) — vectors v and w 



four scalars a, /?, 7 and 5 



Again these shapes remain the same even after Lorentz symmetry violation SO(d —1,1)^ SO{d — 1). 

4.3 Massless gravity 

This is conventional General Relativity with A = B = Q. 

Kinetic matrix (4.1) is actually a symmetric matrix of size '^^''^'"^^ when acting in the space of normal modes: 





00 01 11 


Oa la 


aa ab bb 


00 
01 

11 






1 1 

— \/2^^ll — \/2^^|| 


Oa 
la 








aa 
ab 
bb 


-/cj^ —\/2Ljk\\ —up' 
-fcji -\/2a;fc|| -LiP 







where we denote it by square brackets, keeping ordinary brackets for d x d matrices, like gravity field /i^,^. 
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Its eigenvectors {dx d matrices) and eigenvalues are: 



( 








\ 



















hab 




V 






/ 








... 










a; 






to 






transverse traceless graviton 



transverse (d — 2) — vector 









. . . —OJ 


\ 



























/ 



Stueckelberg (d — 2) — vector 





—w 


k 


... 






/ 





—OJ 


... 




k 














—u 


k 















or 




































v 








) 




I 









temporal or spatial 
Stueckelberg scalar 



W2 





... 





-/,:- 













J 



secondary Stueckelberg scalar 



with 



fc2 


Lok 







\ 


ivk 















A± 














A± 












A± ; 


A± 


= id- 


-3)(- 







two mixed traces, 
spatial and space — time 



= (d - ifoj'^ - 2{d^ - lOd + n)uj'^k'- + {d- Ifk" = {d- - k'y + 16{d - 2)oj'^k 



2^4 



n2/,2 r2N2 



.21:2 



(4.17) 
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In other words, in massless (gauge invariant) case 

Xgr — ^ J 

A+,=u;2^P, (4.18) 

^vec — 

and Qi{\) factorizes, in accordance with (4.15): 

Q4(A) = (a^ + {d- 3)(w' - P)A -{d- 2)(w2 ^ p^2^ ^ ^4_^g) 

so that the two Stueckelberg scalars have A^^^ = and the two trace eigenvalues are 



A± = l\^[d- 3)(-a;2 + p) ± sj{d- 1)2(^4 + ^4) _ 2(^2 _ igrf + 17)^2^2 j (4 20) 

4.4 Mcissive gravity in the rest frame 

Rest frame, where k\\ =0 is distinguished because SO{d — 1) symmetry is fully restored in it and there is no 
reason to distinguish between 1 and other spatial directions. Accordingly the kinetic matrix is: 





00 


Oi 


ii 


ij 


jj 


00 


B-A 





-B 





-B 


Oi 





A 











a 


-B 





B-A 





-u^ + B 













u^-A 





. jj 


-B 





-LO^+B 





B-A 



(4.21) 
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Its eigenvectors {dx d matrices) and eigenvalues are: 





( 











\ 



















transverse traceless graviton 




\ 








hab 


J 







( 










... \ 














1 






— A. 












transverse [d — 2) — vector 







1 













V 








J 






( ° 







1 


... \ 





















A 












Stueckelberg (d — 2) — vector 




1 
























) 








( 


1 


0^ 






A 




1 










Stueckelberg scalar 











J 









/ 











uj^-A 







d-2 







\ 








-1 



secondary Stueckelberg scalar 





( -{d-l)B 













\t + A-B 







I 





\t + A-B 



two mixed traces, spatial and space — time 



(4.22) 



or, if we do not distinguish between the 1 and other directions, 









hij 








) 









-{d- l)B 












+ A-B 



traceless graviton 

Stueckelberg (d — \) — vector 

a mixture of secondary Stueckelberg scalar and space — time trace 

(4.23) 

In other words, one has in the rest frame in the massive gravity case the three possible values of the eigenvalues: 



^iw.c. A 



Aor =0^^ — A 



d — 1 times 
(A±m^ times 



(4.24) 



, _ -{d - 2)uj'^ + dB-2A ± ^{d - 2Y{B - uj-^f + 4(d -Ijg^ 





two traces 
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The masses present in the spectrum can be obtained by solving the equations Ai(a;^ = rr??) = 0. This gives the 

mass of the graviton multiplct = A and the mass of the two mixing scalars = ^ j^^2 f{A-B) ^-l^isition 
= has only one solution). Note that at any values of A and B, one of the two eigenvalues that crosses 
the abscissa axis has the negative slope at the crossing point. This means it is a ghost. One may say that, 
depending on relations between A and B, can be positive or negative, in the latter case the ghost do not 
propagate (and is the tachyon, in fact). However, taking non-zero and large enough k"^, one can get a positive 
solution to the equation = (due to the Lorentz invariance only the invariant combination ui^ — enters 
the equation). 

The only possibility to avoid the ghost is to put A = B (the celebrated Pauli-Fierz case), when the mass of 
the ghost becomes infinite (the equation A^^ = has no solutions at all) and goes away from the spectrum. 
Let us return now to the propagator. The coefficients in the propagator can be rewritten as 



ai = 



a2 = 



a3 



II : = 



as = 



ae 



1 



2to2 
1 



1 - a 
k^ ^ 
1 



1 



k'^ + 



id- 
ai + 



1 



{\-a)B 
A{dB - A) 
1 

^2 



m 
1 

¥ 
1 



2 k^+M"^ 



1 



fc2-|-M2 



{d- 
1 



k^ 

d- 



1 d-2 



l)(d-2) 

1-a f 



A;2 + m2 

2 



+ 



1 



m 
"k^ 



+ 



k'^ + M^ 
{d - 2)M2 + rn2 



(4.25) 



M2(d-1) \k'^ k'^ + M'^ 



Here we again encounter two different dispersion laws: k'^ + w? = and k"^ + M"^ = with 

™2 A 71^2 A{dB — A) 

= ^' ^ = {d-2){A-B) (^-2^) 

There is also a fictitious pole at fc2 = that contributes only a ^ 1 and makes the transverse part of the 
propagator. Therefore, it cancels with the conserved energy-momentum tensor. 

As we already discussed the M"^ mode in some channels (coefficients ai) behaves as a ghost, i.e. enters 
with minus sign as compared to the mode, but this does not happen in the channels attached to conserved 
currents. It, however, enters the radiation, (4.4) in the coefficient a2- Note that this coefficient cancels in the 
other distinguished case when the both masses are equal, 



m 



2 = when A = 2B (4.27) 



The receipt to remove the ghost we mentioned is to bring its mass to infinity. In this Pauli-Fierz limit, the 
coefficients (4.25) becomes 

'I- a 1 

oi = 



Ct2 



as = 02-1- 
ai = 
as = 



2m2 V A;2 w? 

1 1 
{d - l)m2 fc2 + m? 

(1-a) 1 



d-lTO2fc2 

^ ' (4.28) 
11 ^ ' 

"2 ' k'^+m? 

1 



(d- l)(fc2 + m2) 
I d-2 1 l-a/m2 d-2 1 



- 1 fc2 -h m2 V k^ d-lk"^ 
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Interaction of two stress tensors looks in this case especially simple if a = 1: 



1 



^2 _|_ 



1 9 



d-2 



(4.29) 



where k^^T/^v = tv For a conserved stress tensor t = 0. 

However, beyond the quadratic approximation ghosts show up through the Boulware-Deser instability in 
curved backgrounds even in the Pauli-Fierz case (it is enough to consider quadratic perturbations but near the 
non- trivial metric g^^ = prj^^)- 

4.5 Generic frame 
The kinetic matrix is now 





00 


01 


11 


Oa la 


aa 


ab 


bb 


00 


B-A 




-B 




-kf^-B 
— \/2a;fc|| 




-kl-B 
—\piijjk\\ 


01 




A 








11 


-B 




B-A 




-uj^ + B 




-u^ + B 


Oa 








fcj^ + A cuk^i 








la 








uk\\ to^ — A 








aa 




— V^wfcii 


-LO^+B 




B-A 






ab 










u?-k\-A 


bb 




— \/2a;fc|| 


-LO^+B 




-u? + fc2 + B 


B-A 



(4.30) 



Lorentz-invariant eigenmodes are much simpler: they are obtained by Lorentz transformations from those 
in the rest frame. Accordingly the eigenvalues are 



Xgr = iiP" — k^ — A 



d—\ times 
(rf+l)(d-2) 



times 



(4.31) 



± _ (d - 2)P + rfS - 2A ± ^{d - 2)2(fc2 + Bf + 4(d - \)B'^ 
A<,„ — 



two traces 



and the corresponding eigenmodes w'"^ are 











. . . —a; . 
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fc|| 






2 ^^11 






— w 


fcll 







1 







v 
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/ 







d-2 






> 









Kb 









... . 


•• \ 








CO 




-^11 


u 
















/ k^ — 



d-2 



d-2 



kf, -Lj^ 



d-2 I 



d(d-Z) 



d-2 
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/ A±+A + (rf-2)a;2 V2{d-2)wkii 

V2{d-2)uk\\ -A± - A+(d-2)fc|| 



V 



-Xt - A 



-Xt - A 



respectively. As usual, we assume here that the spatial momentum is directed along the first axis. 



4.6 The origin of DVZ discontinuity 

What does this mean? The Newton potential, describing interaction between Too-components of the stress 
tensor can be read off from formula (4.3) and looks like 

p-mr 1 / p-mr -Mr\ 

U{r) --J^- .... - 2)-^ + (4.32) 



with 



Whenever 



r-rf-s (d-l)(d-2) 



-4. S - 0, potential V{r) -^lU - J—\ ='L-1. (4.34) 



which is also the massless- gravity value, except for the special PF case when simultaneously 

^ ~ ^ and M ^ oo : then V{r) ^ ( 1 - = • (4.35) 



This is the DVZ discontinuity. 

In other words, this is the discontinuity due to the different limits of the coefficient in the propagator, 

1 ^' + 3^ 1 ( d-2 



fc2 + A(rf-2)fc2 + ^^^ (d-l)(d-2) U' + m2 ^ fc2 + M2 



(A-B) 

11 A-B 



{d - 1) fc2 
1 1 



0, (4-36) 



(d - 2) fc2 

In formal terms, the story is about the limit of a function 

ax + by 
cx + dy 



A,B^O 



(4.37) 



as x,y —>■ 0: the limit depends on the ratio x/y. It is important that this ambiguity is intimately related to the 
singularity of the function along the line cx + dy = 0. 

In physical terms, in PF case the ghost has infinite mass and decouples, but in generic situation its mass 
tends to zero along with graviton's mass, and they both contribute to the Newton potential. When masses 
exactly coincide (in particular, vanish) ghost simply subtracts/adds to the graviton, but in different way for 
different a/j-structures. In particular it adds in a^, which controls the Newton potential. 
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5 Abandoning Lorentz invariance 



Lorentz violation as a phenomenologically viable possibility was suggested by V.A.Kostelecky and S. Samuel in 
[27] (sec also [28]) already long ago, but it acquired enormous attention quite recently, see [29] for incomplete 
lists of papers about its possible role in particle physics and cosmology. This paper is focused on theoretical 
rather than phenomenological aspects of Lorentz violation, which were also addressed in [29] and partly reviewed 
in [13]. 

5.1 Generalities 

Lorentz violation implies that the global symmetry group SO{l, d— 1) is broken down to SO{d— 1). Surprisingly 
or not, this leads to rather drastic changes in the structure of the theory, revealing all the results of breakdown 
of gauge invariance which remained hidden in Lorentz invariant case. In particular, since different reference 
frames are no longer equivalent, the SO(d — 1) symmetry of the spectrum in the rest frame is broken down to 
helicity symmetry SO{d — 2) in all other frames, what gives rise to highly non-trivial dispersion relations for 
the elementary constituents (different polarizations) of the gravity field. These slightly unusual features remain 
obscured in the case of massive vectors and come out of the shadow only for tensor fields. 

Prom h^v and fc^ one can make: 

• 4 /i-linear SO{d— l)-scalars: /iqo, ha, kihoi, kikjhij, 

• 2 /i-linear SO{d — l)-vectors: hoi, kihij and 

• 1 /i-linear SO{d — l)-tensor: hij 

and thus 14= ^ + 2^ + 1 = 10 + 3 + 1 /i-bilinear structures which contribute into the propagator V^j,u^ai3- The 
current- current interaction then looks like 

" IJ,v,a0-i- — 

= aio{kiToif + aii{kiTijf + a2o7oo(fci%7ij) + a2iTii{kikjTij) + 

+a4oTQi + a4iTy- + feTpo + flBoTooTji + ar,i{Tiiy+ (5.1) 
+ciroo(A;iroi) + C2Tii{kjToj) + C3Toj{kiTij)+ 
+a6o{kiToi){kikjTij) + aei{kikjTij)'^ 

Our notation takes into account that each structure with coefficients Ofc in the Lorentz invariant case is now 

split into two and ako, (iki arc the two independent coefficients. The coefficient h corresponds to what was a 
combination of a4 and as and Cj correspond to new emerged structures. Note that 03 is absent at all here, since 
we no longer keep the parameter a (in the previous section we saw it played no important role in massive cases) 
and define the propagator just by the equation KV = E. This would correspond to 02 = 03. 

Two of these structures, [khkY and [A; ft, fc] [/eft,] do not appear in the Lagrangian, because they are more than 
quadratic in momenta. Other structures enter the Lagrangian with adjusted coeflacients so that the Lorentz 
invariance is not broken in fc^-terms. By now conventional parametrization of quadratic Lagrangian with 
manifest Lorentz and general covariance violation in the mass matrix is 

I (k^kaTjffu + k^k^rjau + k„kaV/3n + k^kpria^ - fc2 (rj^,ar]^i3 + Vi^aVnf3^ - 2 (k^,k„r]ai3 + kakpri^^ + 2k'^r]^^rjai3^ Wh°'^+ 

+2mlhlQ + 4m?ftgi - 2m^ft|j- + 2mlh% - Amlhoohu (5.2) 
Lorentz invariance is restored provided 



ml = 


-B- 


-A, 


mf = 


ml 


= A, 


mj = 


ml 


= B 
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The propagator (5.1) is obtained by inverting the square matrix in the action, 



- 


00 


01 


11 


Oa 


la 


aa 


ab 


&& 


00 


jrin 

U 




—ml 






—ku — mf 

II 4 




—fen — m? 

II 4 


01 




mf 








— V2i-(;fc|l 




— V2wfc|| 

1 1 


11 


-ml 




9 9 






-w^ + ml 




-w^ + m| 


Oa 










a;fc|| 








la 










- ml 








aa 


-fc^ - ml 


—\/2iok\\ 


-w^ + ml 






9 9 

ml — m2 




-w^ + fc^ + m| 


ab 














— A;j^ — m2 




bb 


ml 


— V^w/cii 


— + ml 






— + fcj^ + ml 


9 9 

m^ — m2 



(5.4) 



or, in the rest frame. 





00 




ii 


ij 




00 


ml 










-ml 


Oi 
















ii 


-ml 





2 9 

m3 — m2 





— + m| 


ij 











- m2 





. n 


-ml 





-w^ + m| 





ml — ml 



(5.5) 



For special values of masses the theory acquires some residual gauge invariance [11]: 

— > + C{x^,t) : if mi = m2 = ms = m4 = 

t^t + C°ix\t) : if Too = mi = m4 = (5.6) 

x'^x' + C{t): if mi = 

PF gravity corresponds to mo = 0, mi = m2 = m| = ml. In general the graviton mass is r/72 and the 
corresponding transverse modes are always split from everything else, while the rolo of the other mass- 

parameters is to control non-trivial dispersion relations for the other constituents of the gravity field and their 
severe (and generically inseparable) intermixing. A lot of these complexities disappear if mo = and such model 
exhibits only minor deviations from the intuition developed in experience with Lorentz invariant theories. 

5.2 Rubeikov's approach 

In the pioneering paper [10] V.Rubakov made the first attempt to diagonalize the action (5.2). 

He expressed /i^^ and the action through transverse fields (xij, Wjjj "^i ''') and additional Stueckelberg fields 
(sj, V, a) and, next, through the gauge invariant vector Wi and scalars $ and t: 

hoo = V', 

hoi =Ui + kiV, kiUi = 0, (5.7) 

hij = TSij + Xij + (kiSj + kjSi) + kikjCT, kiXij = kjXij = Oj X^^ ~ 0, fcjSj = 
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Since gauge transformations h^i, — > h^i, + k^'E^ + k^E^ with Sq = ^o, = + hC and ki^i = act on these 
fields as 

Xij ^ ^i^j ~^ ^j^ii 
Ui Ui +U)^i, 
Si ^ -\- ^i, 

V'^V + 2a;Co, (5.8) 
V ^ V + u)C, + (,0, 
CT (7 + 2C, 
T ^ r 

the gauge invariant vector and scalars are'' 

Wi =Ui — LOSi, 

$ = V - 2a;t; + w^tr, (5.9) 

T 

and the action for conserved stress-tensor T^^, k^T^,^ = is 

-k\l -id- l)id - 2)oj\^ + 2k\w'i -id- 2)$r + (^-2X^-3) ^2) ^ ^. .y. . _ 3^.^^. ^ ^^^^ ^ ^y.. 

+mo^2 + 2mi(u- + Pu^) - ma^x?,- + (<^ - 1)^^ + 2fc^(s- + ctt) + k'^a'^^ + mliPa + (d - 1)t)2 - 2m4V'(fc^cr + (d 

(5.10) 

The first line is quadratic part of Einstein-Hilbert action and it contains nothing but gauge invariant fields. 
In general relativity the scalar $ does not have a kinetic term and does not propagate. The second line 
contains gauge and Lorentz-violating terms and depends on all the fields, in particular, provides kinetic terms 
for Stueckelberg fields Si, Vi and a (if TOq and m\ arc both non- vanishing) . 

It is rather straightforward to diagonalize this Lagrangian and analyze its particular eigenvectors: modes 
of the gravity field. This was done in [10] under a strongly simplifying assumption mo = (which also 
guarantees the absence of ghosts), the same kind of analysis in general situation being rather tedious, and only 
the Stueckelberg sector was studied in [11] for mo ^ 0. In what follows we use a slightly different technique, 
which also has an advantage of being easily made algorithmic and thus allows one to make tedious calculations 
in a systematic, way with the help of a computer. In fact, as in the previous section, we use two different ways of 
analyzes: both though manifest constructing the propagator and through immediate diagonalizing the kinetic 
matrix. 



5.3 Pcirticle content of the theory 

Gravitational field is described by symmetric dx d matrix 

As in the Lorentz invariant case, the ''''^'^^^^ components of the gravitational fields split into t'^+^K^' ^) 
components describing a gauge-invariant traceless graviton, one scalar trace and d components of gauging 
(Stueckelberg) vectors, which is in turn split into space-time transverse vector and the secondary Stueckelberg 
scalar, sec cq.(4.9) above. In the rest frame (which exists since the fields are massive) these ('^+iK'^~^) com- 
ponents are all degenerate, being related by the action of SOid — 1) symmetry group. However, if Lorentz 
invariance is broken, degeneration is lifted for fc ^ and the components further split into transverse 

graviton, transverse vector and a scalar, which are representations of the "helicity" group SOid — 2) acting in 
the hypcrplane, orthogonal to k. 

^In order to make a contact of these notations from [10] with notations of the paper [13]: 

T = 2ip Xij = Si = -Fi a = 2E m = Si v = B ^ = 24> 
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These degeneration properties are most concisely reflected in determinant formula for the kinetic matrix 

where ^ 

-D^ec = —m\kQ + m|fc^ + m\m\ = m\ { — Aiq + ^^fc^ + 7712), (5-12) 

Dgr = -kl + P +ml (5.13) 
and £) is a sophisticated expression of power 4 in fco and |fc|: 

D = {d- 2)mlmlk^ +{{2d- 4)m^m^ - {2d - 4)m^mi + {2d - A)m\ml - {2d - 4)m^) klP+ 

+ {{d - 2)mlml - {d - 2)m\m'^ 
+ {{d - l)mlml -{d- 3)m|m§m? - {d - l)rr7|m^m?) fco+ (^-1"^) 
+ — 3)mom2mi — (d — 3)7713 TOqTOi ~ (2c? — ^)rn\rn\m\ + (d — 3)to4TOi^P+ 
+(d — l)m|mom3mi — m^mlmQ — {d— l)m\m\m\ 

For d = 4 (5.14) gives 

£> = 2mlm\kQ + {-Amlml + Amlml + Arn\rnl - Aml)klP + {2mlml - 2m\ml)P+ 
+{—rn^rn^rn{ — Sm^mQmi + Sm^niiJkQ + {rn^m^ml — rn^m^ml — ^w^m^mi + m\m\)l?+ (5.15) 

,o2222 422 n224 

+6m,2mQm^mi — m2mimQ — 3m2mim^ 
At k = our -D^ec turns into Dgr and D decomposes into two factors, one of which is also D2: 

_ fc=0 2 r-) 

D H DgrD,r, (5-16) 

Dtr = ml {mlml + {d— l)m\ + (d — 2)mQa;^ — (d — l)mQm3) 

Clearly a drastic simplification occurs also in the case of mo = [10], when D turns into 

Dmo=o = ^m\{m\ - ml)kQp + 2ml{ml - ml)P + 'imlm\kl + m\m\{m\ - 4m|)P - 2,m\mlm\ (5.17) 

One more implication of (5.11) is that only degrees of freedom from the first braces in expansion (4.9) are 
dynamical fields, even after violation of the Lorentz and gauge symmetries. 

5.4 Normal modes for gravity fields 

We manifestly describe the normal modes only in the rest frame in this case, since, formulas in the moving 
frames become very involved and non transparent. Note that Lorentz transformation can no longer be used to 
deduce formulas in the moving frame. Therefore, we explicitly present and discuss, at least, eigenvalues in the 
generic frame. Note that the rest frame analysis is already enough to see ghosts. 

Rest frame, where fc|| = is distinguished because SO{d — 1) symmetry is fully restored in it and there is 
no reason to distinguish between 1 and other spatial directions. 
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Its eigenvectors {dx d matrices) and eigenvalues are: 
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\ 




- ml 















transverse traceless graviton 




\ 








hab 


J 





/ 




Lu — m% 



1 

V ••• 

/ 





transverse (d — 2) — vector 



m\ 



1 



Stueckelberg (d — 2) — vector 





1 





1 





m\ 




1 










\ 












Stueckelberg scalar 





( 











- ml 







d-2 







\ 








-1 



secondary Stueckelberg scalar 





' -{d- l)ml 







A± - ml 







\ ^ 





two mixed traces, spatial and space — time 



or, if we do not distinguish between the 1 and other space directions, 















mi 














-{d- l)ml 








A^c - ml , 



(5.18) 



traceless graviton 
Stueckelberg (d — 1) — vector 
a mixture of secondary Stueckelberg scalar and space — time trace 



(5.19) 

In other words, in the massive gravity one has the three possible kinds of the eigenvalues in the rest frame: 

(d+l)(d-2) 



Xgj- = Up' — m2 



m^ 



times 



Ati„ 

_ {d- l)ml ml + m2 {d- 2)uj'^ 



d — 1 times 
two traces 



(5.20) 
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The last two eigenvalues, A^^ describe only one propagating mass. Indeed, the zeroes of the two equations 
A^^ = are encoded in the equation A+ A"^ = 2^mQ + — (d — l)m| = 0, which is linear in lo"^. This mass is 
manifestly given by 



Therefore, in this case of broken Lorentz symmetry, similarly to the Lorentz invariant case, there are two 
propagating masses, and = ml, although the dispersion laws are far more tricky in the present case, 
as we shall see in s.5.6 (e.g. some of these modes propagate with another speed of light). Note that the two 
masses coincide, = m^ provided 

[ml — m\)m1 = m\ (5.22) 

Now, similarly to the Lorentz invariant case, one immediately observes a ghost: since the derivative ^^ji _^ 

is negative at that single point where A(w^) crosses the abscissa axis. This ghost may not propagate if the 
parameters are chosen so that < 0. However, then one has to look at non-zero spatial momentum to see if 
this crossing point A^ becomes positive. 

As before, one may try to remove this ghost from the spectrum bringing its mass to infinity, which is 
equivalent to putting mo = provided the both mo and are finite, or either mo or m4 go to infinity. 
Therefore, mo = is, at least, a sufficient condition for ghost free space-time, [10]! 

5.5 Characteristic polynomials 

Now let us analyze the modes in general frame looking at the eigenvalues. Here we present the analysis in terms 
of the Euclidean normal modes (using the Lorentz normal modes this time does not lead to any simplifications, 
since the Lorentz symmetry is violated anyway). 

One of the eigenvalues corresponds to the graviton propagator and is 



Dgr = -\gr = -oJ +ml = k^ + ml, (5.23) 



It describes the '^^'^^ graviton modes. 

The Stueckelberg and transverse vectors are described now by the 2 x (rf — 2) eigenvalues'^ 

A^ec = \[^'' + k'' + ml-ml± ^ (a;2 + p)2 + 2{-uj^ + k^){ml + ml) + {m\ + mlf^ = p2 ± ^P^, 

(5.25) 

which combine into the polynomial 

Dyec = -KecKec = Pi - pI = -mlu!'^ + mil? + mlml (5.26) 
and are described by a characteristic polynomial 

P2(A) = A2-2p2A-£>„ec (5.27) 



•^One can compare these formulas with those for the Lorentz normal modes, 

^ _ fe2 _^ mf + m| ± ^J{ml - ml)'^ + 2{ml - m|)(a;2 + + 

^vec ^ 

One can see there axe no any simplifications in this case. 



(5.24) 
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(5.28) 



The remaining four scalars combines into quite involved quartic characteristic polynomial of the form: 

Q^{X) ^ X* + X^(^{d-3){ui^ - P - ml) -ml-ml + 2m\ - 2m^) + 
+A2| - {^d-2){u? + k^f +u?i^m\ + (d - 3)(-mg - + + m^)) + 
- 2m| + (rf - 3)(mo + - + - 2m|)^ + 

I // 2 2 1 4 o 2 2 o 2 2\ I / J -i \ / 4 2 2, 2 2, 2 2\\ 1 

+ 1 (rnQTO]^ + 7712 — lm,Qm2 — lm^'m2) + (a — 7724 — ra^ra.^ + mj^mg + m^mQ) I > 
+A|(d - 2)(^a;^(mo +mi) + 2uj'^P{ml -ml+ml + ml) + P{ml - ml+ml)^ + 
— l)(m4 — (toq + mf)m|) + {d — 3)(— mfm2 + m^m^l — mQm%)^ + 
+P{2{d — 2){m\m\ — m^m^) + {d— 3)(mQm2 + m^m"^ — m^ml — m\m1 — mlm^ + m'^ + 

1/0 222 24 24, /J ^\^ 111, 42 24 222, 22 , 

+ 1 2mQmim2 — TnQm2 — mim2 + (a — 1 j {mQm2m^ + m^mi — m2m4 — m^mim^ + mim2m^) 1 > + 
+ ^{d — 2)^ — uj^mQm\ + 2uj'^l?{m\ — m^m^ + m^m^ — m^m|) + k^{m,1m,2 — mfm|)^ + 

(j^d — l)(m^^m\ml — m\m\) + (d — 'i)m^m\m\^ + 

+J?{2{d — 2)m\m\m\ + {d — 3)(— rnQm^ml + mQm\m\ — m\m\)^ + 

+ {m\m^rn[ + {d— l){m\m\m!{ — m\m\mQm\)^ | 

so that Qa{0) = —D. This quantity is definitely the same in the Lorentz and Euclidean modes, since this comes 
from the determinant of the kinetic operator in the scalar mode sector (similarly, are invariant Dgj. and Dyec)- 
In the rest frame, where = 0, we definitely returns to formulas of the previous subsection: 

\ 2 2 

Xgr = OJ — m2, 



A+ , = uj^-ml (5.29) 

and (34(A) factorizes, in accordance with (4.13): 

Q4(A) = (A-m?)(A-(a;2-m^)). 
•|a2 - x(ml -ml + {d- l)ml -{d- 2)uj^^ - (d - 2)Lu^ml - mlmj - (d - 1)to| + {d- l)mgm^| 



(5.30) 



Therefore, one can see in the moving frame that, of the ^) modes in the rest frame with A = — m^, 

remains graviton modes, while d — 2 modes along with d — 2 Stueckelberg modes (with A = ml) compose 
the -D^ec) and one mode with X = co^ — m\, one mode with A = m{ and two trace modes altogether compose D 
(or Q4(A)). 

5.6 Ghosts, tachyons and others 

With the characteristic polynomial in hands, we can discuss the properties - and peculiarities - of the spectrum 
of Lorentz-violating gravity. Basically, there are three different important peculiarities: when there is a ghost, 
when there is a tachyon and when the speed of propagation of a mode differs from the light speed. First of all, 
the tensor sector is healthy of all these phenomena. 

Superluminal propagation. One can most directly observe the superluminals in the vector sector of the 

theory. Indeed, the dispersion law D^^c = implies that the vectors propagate with the speed ^ia This speed 
can be larger or less than the speed of graviton (which is 1) depending on the relation between masses. 
However, the main peculiarities are related to the scalar sector of the theory. 
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Ghost. We already discussed the appearance of a ghost in the scalar sector in s.5.4 and concluded that it may 
have infinite mass (i.e. disappears from the spectrum of the linearized gravity, but can easily come back beyond 
quadratic approximation as a Boulware-Deser mode), provided mg = 0. Otherwise, there may be a ghost, at 
least, at = and > 0. In other words, there will be a mode which is constant in space and would grow 
in time (simultaneously in the whole space). Still, this is not a ghost propagating in space. 

In order to move slightlj' away from the P = point, i.e. consider propagation in space, one can find a 
solution to (5.28) at small P. The limit of P is smooth, and one would come to a ghost again unless mi = 0. 
Indeed, one expects a singularity in this latter case, since the determinant (5.11) is proportional to mi at zero 
momentum. 

Let us make a closer inspection of this particular case. In this case the dispersion law looks quite strange 

D ~ oj'^P = (5.31) 

This corresponds to infinite speed of light for the vector mode and leads to a very exotic excitation, being a 
carrier of instantaneous interaction = 0. One can actually describe the slightly-virtual (small values of w^) 
instantaneon analytically: the corresponding eigenvalue is 

. 2(d-2)A.a;^P 

(d-2)(m§-m2)fc4-[(d-3)A + 2(rf-2)m2m2]fc2 + (d-l)m2(m2m§-m|)-m2m4 ^ ' 

(5.32) 

where 

A = mliml - mi) - m| = (5.33) 

This excitation is simultaneously a ghost, dX/du)^ < on shell, when P lies in between the zeroes of the 
denominator in (5.32). 

Thus, in the specific case of mi = the ghost is related with excitations with the specific dispersion law, 
and is not a particle-like ghost. This is probably the reason why it was not recognized as a ghost in [11]. It 
would be interesting to better understand physical implications of this instantaneon. 

Tachyons. Now we return to our discussions of the ghost-free regime at mo = 0. In order to see if there are 
tachyons, one has to put w = and look for real solutions X{k) = 0. Note that this is not the same as to look 
for a mode with negative mass square because of the tricky dispersion law. Indeed, the mass describes the pole 
of the propagator in oj at zero k, while the tachyon has to do with its /c-dependence. 

In order to guarantee the absence of tachyons in the vector and tensor sectors, one has to require mf > 
and 1712 — 0- Then, it is again enough to look at the product of all 4 scalar eigenvalues, which is D. Thus, the 
tachyon is absent as soon as there is no real-valued solution of the equation 

-£)(a;^P) =pfc* + r?P + C = (5.34) 

a)=0 

where 

p= {d — 2)(mim2 — mlml) 
rj = 2{d — 2)m{rn^ml + (d — 3)(— mgmim^ + m'^m\m\ — m\m\) (5.35) 

> 422i/j i\/242 222 2\ /j o^ 2227l/r2 

Q = mjmgm]^ + (a — lj{m2m^mi — in2m^inQmi) = —{d — Z)rnQm^m2M 
The tachyon is absent either if the discriminant of (5.34) is negative, 

rf - 4pC < (5.36) 
or if both solutions \ (^—rj ± \J'if' — 4pr?^ of the quadratic equation, (5.34) are negative 

7? > 0, pC > (5.37) 

If ncithca- of these conditions is satisfied, there is a tachyon in the spectrum. In order to have a theory both 
without the ghost and the tachyon, one can put m^ = and require that 

m\ > ml, 2{d - 2)m\ > {d- 3)ml, m| > (5.38) 
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(these are the conditions obtained in [10] in = 4) or 

[2{d - 2)m^ -{d- 3)to^] ^ < 4((i - l){d - 2) 'm\{'m\ — mg) (5.39) 

We assumed here that mf , and m\ are non-zero. 

One can also consider the border cases. If one of the masses mf , rn| or m\ is zero, C = which means that 
there is a massless mode in the spectrum. Now one has to differ between different cases. 

If ml = 0, there is a tachyon unless also = 0. In this latter case, the speed of light of the vector mode 
becomes zero, and the dispersion law 

£) ~ fc* = (5.40) 

implies the mode does not propagate in time at all. 

There is also a possibility of m| = that leads to a non-propagating mode as well, with the same dispersion 
law (5.40). 

The last border case to consider is p = 0, i.e. = m%. Then, the tachyon is absent if > 0. This means 

2(d- 2)m| > (d- 3)mi > (5.41) 
In particular, if the equality is realized in these formulas, 77 = and 

2{d-2)ml = {d-2,)ml (5.42) 

the dispersion law acquires the form 

2(d - 2){ml - ml)Lo^P -{d- \)m\m\i^'^ + {d- \)m\mlm\ = (5.43) 



5.7 Dispersion relations 

One of spectacular puzzles of massive gravity is emergency of non-trivial dispersion relations a; = e(|fc|). This 
looks puzzling because usually does not happen in the theory of some N scalar massless fields perturbed by an 
arbitrary mass matrix: 

N N 
a 

)+ Yl ^o-b(l>a(l>b (5.44) 
a=l 0,6=1 

In general kinetic and mass matrices define two quadratic forms which can be simultaneously diagonalized, but 

diagonalization of Mab breaks down diagonal form of the field-current coupling. This phenomenon is well known 
as the Kobayashi-Maskawa mixing in the Standard Model of elementary particles [30] . Characteristic equation 
defining the eigenvalues of such quadratic Lagrangian without currents. 

Dm = det (k^dab + M^b) = (5.45) 

is actually a product 

N N 

Dm = J{ {k^ + ml{M)) =l[(^-u;' + el{\k\)) (5.46) 

a=l a=l 

so that dispersion law is the standard relativistic uj = ea(|fc|) = \Jk'^ + m^. 

If Lorentz symmetry SO{d —1,1) is broken down to SO{d — 1), there are in general three different matrices: 

N N 
a=l a, 6=1 

and characteristic equation, defining ci;(|A;|), is more sophisticated: 

Dn,m = det f - u^dab + Nabk'^ + Mab) = (5.48) 

NxN \ J ^ ' 

While Dn^m is still a product like the last formula in (5.46), the roots ea(|^|) can now be highly non-trivial 
functions of the space momentum k. This would explain the origin of non- trivial dispersion relations, but the 
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problem is that in massive gravity one docs not introduce any non-trivial matrix Nab 7^ Sab- all Lorentz violation 
is concentrated in massive matrix and docs not affect the kinetic term! This seems to imply that nothing more 
than a KM mixing can occur with no severe change to dispersion relations - but this is actually not the case, 
as we see in eqs.(5.12)-(5.14) and in s.5.6. 

Resolution of the puzzle is in the concept of Stueckelberg fields: when mass matrix breaks some gauge 
symmetry, it gives rise to kinetic terms for the newly revived gauge degrees of freedom. From the point of view 
of above scalar theory this looks strange: if massive matrix involves more fields than the kinetic matrix, 

N N+n 

(k^cg + Ja(t>a) + ^ Mab(i>a4>b (5.49) 
a=l a, 6=1 

then the extra n fields 4>a should be represented as derivatives, (pa' = J2b'=i ^■t'^a'w'^b'- Then the new kinetic 
term involves new Stueckelberg fields ipa' ■, breaks Lorentz invariance and we come back to the situation described 
in (5.48), where non-trivial dispersion relations are of no surprise. The problem is that above substitution (j> ^ if 
looks ad hoc and it is actually justified only when (pa' describe pure gauge degrees of freedom. 



5.8 Manifest expression for the propagator 

Now wc come the current-current interaction. To this end, we need to know the propagator. Coefficients in the 
propagator in 4 dimensions are: 

aio = ^ ((4^4 + 2rn\mQ — ArnQm\)u)'^ — im'^mlu'^k'^ + (2m| — 2rn\'m\)k'^+ 

-|-(7m2mgm3 — 7m\m\ — irn^m\)Lij^ + (4m2m4 — Tn^m^ — m^m^ + vi\m^ni\)k'^ + {Sm'^m'^ — 3mQm1m2 + ni^n 

~2{-uj^ + P + mf) -LO^ + k'^ + mj 

'^ii = =: ZT^ = ~2- 



m\{—uj'^ + k'^ + ml){—uj^ + ^k^ + m^) Dy^cDgr 

{2{ml — 2m|)a;^ -|- m\m\ + mlml — 3mf mg) 
a2o = -2 — , 

2 f / 2 2 , 22io22 o22 o4\2i/-22i 22 2 2\r2i 

a2i = — — — I (m^nii + rnQm-^ + Irn.^rn^., — Zm,2m,^^ — Zm^juj + (m^m-^ + rn^m^ — TOj^TOj)/? -|- 

J-^ gf -/-/ 

i/ri24i 222 222 o22 2\\ 

+ {2mim^ + m2'm-^^mQ — 'm^'m-^^m2 — 2rn^mQrn-^) 1 , 

aeo = — — ((2m2mQ — Arn^rn^ — 2m\m\ + Am^oJ^ -\- {2'm\m\ — 2m\m\ — 2rn\rn'\)l?' + 

j-J 1^ (-J LJ \ 

, I 222, 222 42io222io222 24 o2 4\\ 

+(~"^2"^i"^o + TOsTOq™! ~ ^2*^0 + 3m2momg + 2m^mym2 — rriim^ — 3m2m^) I , 
aei = — — i{^rn\mQ + 2m\ + ?,m\ - d,m\m\ - ?,mQm\)uj^+ 

^vec^ gr ^ ^ 

+ {?>m\mi + imlml - 8rn| - Amlm\ - 2m\ - Am\m\ + 2m\m\)uj'^P + {Am\m\ - 2m\m\)k'^+ 

-|-(8mim|m| — Am\rn^m\ — Am\m\ + Am\mQ'm\+ 
+&rn^rn^rnl — 2m'^m\ — dm^m'l — rn^m\ + 2m\m\ — Smfrn^ + mfmQ)Lj'^+ 

i/'r>24,o222i 222 o222 o222 42,j42 02 4x7*2 1 

+(2m-im^ + 3m2m^mi + mQm2mi — 2m^mQmi — 2'm^m2mi — m2mi + AmiUi^ — 2m2mi)k + 

1/1 O/IO 010 010 /lOO A A \ \ 

+ {2rnim^ — 2m2TOx™4 + TO2"^i"^o + 'S'm2mi'm^ — 2mim^mQ — ■m2m-^) ] , 

{-2m\)ijj^ + 4(m| - ml)oj'^k'^ + {m\m\ + 3m\rn{)uj'^ + {mlml - m\m%)k'^ + {m\m\ - 3m\m\ml) 

D 

-up' + m\ 
040 = 2- 



D, 

041 = - 



vec 
1 

(5.50) 
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051 = ^ ({mQm1)uJ^ + {2rn\'m\ + 2m^ml — 2m\ — 2m^m\)LJ^k'^ + m\{m\ — ml)k'^+ 

J-J gy i-J ^ 

2m\ ^ - rn\uP' — (mg - ml)P + m'^'m?^ 



050 = 



ci = 



D 



C2 = 4a; 



D 

(ml + rri^mQ — m§TOo) 



C3 



D 

4a; 



■^vec 

(5.51) 

Restoring the d-dependence, one can see that only a few coefficients slightly depend on the space-time dimension: 

I ((2d - 4)to| + id- 2)mlml - {2d - A)mlmfj u}^ - {2d - A)mlmluj'^P +[{d- 2)m\ - {d - 2)mlmfj k'^+ 
+ (j(3d — 5)m2mQm3 — (3d — b)m\m\ — {d— l)mQm|^a;^-|- 
+ (j^2d - A)m\ml -{d- 3)miml -{d- 3)mlmi + {d- 3)mlmlml^P+ 

+ (j^d — 1)to2TO4 — (d — l)m^m\m2 + m^rn^ |, 

-2{-uj'^ + P + mj) -u!^ + P + ml 

'^ii = ~ z::2-^ = ~2 



020 



m1{—L0^ + fc^ + m2)( — + ^fc^ + m^) DyecDgr 

{{d - 2){ml - 2ml)u'^ + {d- 3)m?ml + mfmi - (d - l)mfm|) 



D 

2 



021 = j: — p: ({mlml + mQ-mf + 2m\m^ — 21X1^1X1^ — 2m'^u)'^ + {rn\m\ + m\m\ — m\m2)l?+ 

J-J gy -lV V 

+{2m\m\ + m%m\m'^ — rn'lm\m\ — 2m\mQrn\')^ , 

O60 = -j: ({rn^rriQ — 2m\m^ — m^ml + 2TO|)(d — 2)a;'^ + {m\m^ — m\m\ — m\m\){d — 2)fc^+ 

+(— (d — 3)m2m1mQ + (d — 3)m^mQmi — m^mQ + (d — Vjui^rn^rn^ + (d — 2)m\m\m2 — (d — 3)m{m\ — (d — l)m\m'^ 
Oei = — — — — f (2m2mg + + Am\ — ^m^m^ — Arn^rn^){d — 2)a;'^+ 

^vec^ gr ^ ^ 

+{imlml + Amlml - Am\ - 2mlml -m\- 2mlml + m\ml){d - 2)a;^fc^ + {2mlml - mlml){d - 2)fc^+ 

-|-((4d - 8)mlmlml - (4d - 12)mlmlml - (4d - 12)m|mf + (4d - 12)mlmlm\+ 

+ {2d - 2)mlmlml - 2mlml - {2d - 2)mlmi - (d - 3)m2mi + (2d - 6)miml - (d - l)miml + (d - 3)mtmoV^+ 

-|-((2d — 6)mlm\ + (d — \)rn\m\rn\ + (d — 3)m^rn\m\ — {2d — 6)m?^rnlmi — {2d — (!>)m\m^m\ — m2ra\+ 

+ (2d - A)miml - (d - 2)mlm\)1?+ 

-|-((2d — &)m\m\ — {2d — Q)m\m\m\ + (d — 3)m|mimQ + (d — l)m\m\m\ — {2d — Q)m\m\mQ — mgm^J)^ , 

(5.52) 
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+ {{d — 3)m|mi + {d— l)mlml)u!^ + (d — 3)(m3mi — m\m%)k^ + {m2m\ — {d— l)m\rn%m\) 

040 = ^- 



Dyec 
1 

an = - — , 

_l_y Q rp J_y \ 

,/42 222 222x21/222 o222i 42 222\r*2,/2222 22 

+(m4mi — m2TnQmi — m^mQiniju + (mQm2mi — 2m^m2mi + m^rrii — m^mQmijk + {1712171^17131711 — m2mi 



050 



Cl 



D 

-A{d - 2)w{m\uP' + (mg - m^fc^ - mlmf) 



C2 = 4a; 



D 

{m\ + m^mQ — mlmg) 



C3 



4a; 



(5.53) 



5.9 The case of conserved currents: 

For conserved stress-tensor k^T^"' = ojT^" + hT^^" = T^'^k^ = (5.1) turns into: 

Vi,^,a0T^''T°''^ = (aiow^ + a2ow^ + 6 - cia; - oeow^ + a6ia;*)To^o+ 
+{anu)'^ + - C3UJ)TQ^ + (a2ia;^ + 050 - C20j)TooTii + aiiT^^ + 051 (Tjj)^ 

and in d = 4 

1 

tensor propagator : = —-p-, 

vector propagator : [auLO + 040 — C3L0) = — I 4a; Dgr — Dy^c + {rrii — 777.2)7712 ) 



(5.54) 



(5.55) 



No any essential cancelations happen in these expressions, and answers for the scalar channels look very lengthy 
and involved. 

5.10 No DVZ jump 

According to our treatment of DVZ discontinuity it occurs if one removing a ghost from the spectrum via 

bringing its mass to infinity, simultaneously removes its contribution from the static potential, the quantity 
controlled by fc^-dependence. In the Lorentz-invariant case these two things inevitably happens together. On 
the contrary, in the non Lorentz invariant case, the mass of the ghost (=omeg'a^-behaviour) and the static 
potential are unrelated and, therefore, the DVZ jump does not happen. Indeed, if one sends ?Tio to zero, it 
makes the ghost mass infinite and removes it from spectrum via canceling the coefficient in front of a;** in D, 
(5.14). At the same time, the coefficient in front of k^ in D becomes non-zero in this case, therefore, not 
changing asymptotics of the static potential. 
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6 Mixing with extra fields and Kaluza-Klein theory 



Somewhat amusingly, the mixing of gravity with additional field was considered already in the seminal paper 
[2], however there it was used just as a technical trick. Recently this kind of modification of massive gravity 
was re-introduced by S.Dubovsky [11]. Since then the subject attracts a rapidly increasing attention. The idea 
is to add terms like 



and their Lorentz- violating analogues to the quadratic Lagrangian, thus introducing mixing of gravity field with 
something else, what is denoted by tt in these formulas. These 7r-fields can be considered as shifted vector 
or scalar fields (say, Goldstonc fields describing fluctuations near the vacuum expectation values, which cause 
spontaneous violation of gauge and Lorentz invariances). There is already convincing evidence that such mixing 
can substantially soften the strange properties of massive gravity and provide healthy perturbatively-reliable 
models with massive graviton and Lorentz violations. 

Of course, this conclusion is of no surprise, because such healthy theory is well known for decades: this is 
nothing but the ordinary Kaluza-Klein gravity. 

6.1 Example of Kaluza-Klein graviton, d + 1 — 5 '^""^^'^f};!^'^^^'^^ d — 4: 

Kaluza-Klein (KK) gravity is ordinary general relativity in higher d + m dimensional space-time, compactified 
back into d dimensions. From d-dimensional point of view the theory looks as an infinite KK tower of fields 
with different masses, all interacting among themselves. However, in quadratic approximation the fields from 
different KK sectors (with different masses) do not interact and even mix. and one can safely consider each 
sector separately. In this sector we have a massive d-dimensional graviton, which should be completely free 
of any kinds of problems, even gauge invariance (under general coordinate transformations in d-dimensions) 
is preserved. The question is how this can be consistent with seemingly imavoidable pathologies of massive 
gravity, discussed in the previous sections. The answer is that this massive graviton is being mixed with the 
other KK fields of the same sector (with the same masses), and this is the simplest possible argument that 
addition of extra fields can cure massive gravity from all its potential problems. 

Our task now is to analyze massive KK gravity (i.e. a given mass level of KK tower) in some detail in order 
to see how it works. We restrict consideration to one extra dimension, compactified on a circle of radius Rd, 
and express everything in units of Rd- Since even quadratic action for massive KK fields is not widely known, 
we begin with its detailed derivation. 

6.1.1 Quadratic part of KK action in a given sector 

Denote different components oi d + 1-dimensional graviton through h^i,, = /i^^ and (f> = hdd, where = 
0, 1, . . . , d — 1. The discrete momentum in compactified direction is kd = n, n is an integer multiple of inverse 
radius of Rd- We consider a particular sector of a given n, it is not mixed with other sectors in quadratic 
approximation. 

Einstein-Hilbert action in this approximation is 



h'^'^kfji'Kv -l- TT — squared terms 



(6.1) 



or 



TT + TT — squared terms 



(6.2) 



2(fc/i)if - k^hliN - ^{khk)h + k^h^ — > 
2(^^/1"^ + nA^f + 2{k^A'' + ncj,f - {k^ + n^){hl, + 2^ + <^2)_ 
-2{k^,k^h'"' + 2nk^,A'' + n'^(p){h + (p) + (fc^ + n^){h + = 

= {2{kh)l - k^hl^ - 2{khk)h + + n^{h^ - hl^)+ 



(6.3) 



+2 (fc^fc, - k\,) A'' A'' + An {k^h^"" - k^h) A,+ 



+2 {k'^h - {khk)) cf) 




(6.4) 
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which is taken into account by the standard parametrization of Kaluza-Klein (KK) metric, 



e 









: 






\ 





\ 



A, 



d-2 



with (j) = g^f (p. After this shift we get the quadratic part of KK action in the form: 

{2{kh)l - {k^ + n^)hl^ - 2{khk)h + (fc^ + n^)/i^}+ 

+2 (fc^fc, - ^Tjf,,) A'^A'' + An ikf.h'"' - k"h) A,+ 



non — hnear terms 
(6.5) 



(6.6) 



+ 



d- If 
d^\~' 



k^^ + 2n(j){2k^Ai^ - nh) + 



d-2 



Now we have a famiUar action in massless sector (n = 0), which describes gravity plus photodynamics plus 
additional neutral Brans-Dickc scalar. However, for n ^ Q the action is still strange: graviton has mass n but 
there is no mass term for the photon and scalar has mass, different from n (worse than that, the scalar "mass 
term" has a wrong sign!). Instead there is severe mixing between all the three fields: h^v, A^ and (j). To 
highlight the problem we can rewrite the last line in (6.6) as 



1_1 ( _ (fc2 + ^2)^2 ^ 2ncl){2k^A'^ - nh)) + 2 



ncj) 



(6.7) 



and especially strange is the last term. 

Of course, diagonalization of this action is not a big problem - and in fact a literal repetition of that for 
massless gravity, only in d + 1 space-time dimensions. Then, one certainly obtains (d-l- l)(d — 2)/2 propagating 
modes with the mass n (which form the tensor multiplet of (d+ l)-diuK;nsional gravity with the d-th component 
of spatial momentum equal to n), d+1 zero modes (corresponding to the Stueckelberg fields) and d{d + l)/2 
non-propagating modes (which involve longitudinal graviton). This result is guessed without any calculations, 
after some prejudices are thrown away. Still, before we proceed to the answer, it is instructive to analyze 
immediate peculiarities of KK gravity. 

6.1.2 Properties of massive KK graviton 

• The last term in the first line of the rf-dimensional action (6.3) implies that the KK graviton corre- 
sponds to the PF choice A = B. 

• However, gauge invariance is not broken, because of the h — A — (j) mixing. Indeed, one can easily 
check that (6.3) is invariant under 

^h^i, = k^S^u + k^^fj^, 



6Ai^ = + k^C, 
S(t> = 2nC 

• Furthermore, there is no ghost, because this is PF gravity. 

• As already mentioned, there should not be any BD instability. 



(6.8 



6.1.3 Diagonalizing KK massive sector 

We can now return to the problem of analyzing the KK Lagrangian (6.6). In fact, it is better to return one 
step back - to (6.3). The peculiarity of KK theory is that the massless (n = 0) and massive {n ^ 0) should be 
handled in two very different ways. While in the massless sector one makes the celebrated shift (6.4), which 
leads to complete separation of massless graviton, vector and scalar fields in kinetic matrix, in massive sector 
one should make an absolutely different shift: 

h^i/ ^ h^i, -\- {k^Ai, + kifAfj^ j^/^^^^ (^*^) 
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Prom the very beginning note the two things: first, coefficients have n in the denominator, thus this shift could 
not be done in the massless sector, and second, after this shift the new field /i^^ is not affected by the gauge 
transformations (6.8) at all - it should be itself considered as gauge invariant. The result of (6.9) on (6.3) is 
spectacular: the last two lines are fully eliminated, i.e. this shift reduces the sector n 7^ to just a single PF 
massive graviton, 

2{kh)\i - Phltj^ - 2{khk)h + fc2/i2 — , 
^-^ {2ikh)l - kX. - mk)h + k'h'} + n'{h' - h%)+ 
+2 (fc^fc^ - k'^ri^^) Ai^A" + An {k^h'"' - k^h) A^ + 2 {k'^h - (khk)) cj) 
n = / (6.4) (6.9) \ n ^ 



[2{kh)l - k^hl^ - 2{khk)h + {2ikh)l - ehl, - 2{khk)h + fc2/i2| + „2^/j2 _ ;^2^) 

This is of course what one could expect from the very beginning: the ^J±m^ de grees of freedom of a + 1- 
dimensional massless graviton can turn either into the ('^+iK'^~2) _ d(d-a,) _|_ _ 2) _|_ x rnodes of massless 
graviton + massless vector + scalar in d dimensions in the n = sector or into the ('^+^)^('^~2) degrees of 
freedom of a d-dimensional massive graviton in the n ^ sectors. There is simply no room for anything but 
a massive graviton in n 7^ sector, thus nothing like massive vector or scalar can exist there in addition to 
massive graviton. 

It is also instructive to look at the same counting from the point of view of the eigenvalues. The field Hmn 
in d + 1 dimensions had a kinetic matrix with ('^+^K'^+2) eigenvalues A, of which exactly d + 1 were vanishing. 
This left ('^+iK'^+2) — (^d+ 1) = '^^'^^^'> eigenvalues - as needed for kinetic matrix of a li-dimensional symmetric 
matrix. We saw in s.4 that all these eigenvalues are indeed non-vanishing - though not all associated with 
propagating particles. 

6.1.4 Compactification of a vector field 

To get a better illustration of what happened in the previous subsection, one can repeat the same trick for a pho- 
ton field: consider a massless {d+ l)-dimensional photon Am and look what happens to it after compactification 
to d dimensions, where it turns into a d-component vector and a scalar (/). The Lagrangian 

{kMkN - (fc^ + n^)5MN)A'^A^ = (k^k^ - {k^ + n^)) A^A"" + 2n(fc,,^^)0 - fc^^^ (g_;LO) 

In the simplest case of d = 1 one has in (1 -|- 1) = 2 dimensions just {kcj) — nAf with the kinetic matrix 

(-'n? kn \ 
kn ) 

and eigenvalues A = and X = k'^ +n'^. Prom the one-dimensional point of view, one has instead a single mode 
{k(j) — nA) with the eigenvalue A = 1. The difference is dictated by normalization: the usual fact for quadratic 
forms. 

A similar phenomenon occurs for gravitons: non- vanishing eigenvalues in d -|- 1 and d dimensions are in 
one-to-one correspondence, but do not literally coincide because of different normalizations. 

6.2 DVZ discontinuity 

Since KK graviton is the PF one, it is not a too big surprise that DVZ jump occurs when KK radius tends 
to infinity and a graviton mode with a given n ^ becomes massless. However, in KK case the discontinuity 
has a simple explanation: there are additional fields, A^ and (f) and they also contribute to the interaction of 
the stress tensors. Discontinuity is exactly the contribution of these extra fields. 

The Born interaction between two stress-tensors through exchange of KK graviton from the given-n sector 
can be immediately read out from the massless case in (4.36) by making two changes: d ^ d+1 and ka — > n. 
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Wo consider only the interaction between conserved d-dimensional stress tensors, while T^a = Tm = (actually 
this does not affect the formulas too much). The result is: 



and in the massless limit n — > (i.e. Rd — > oo) we obtain 

massless limit of KK gravity mode : ^ i^^v ~ ^ ^ \ ^^^ i^-^"^) 

what is different from the answer for ordinary massless gravity in d dimensions, 

massless gravity : ^ (t^, - ^^T^) (6-13) 

As already said, this is not a big surprise, because to KK answer the other fields are contributing. Moreover, 
(6.12) coincides with the contribution of the massless KK sector with n = 0, and there the contribution of the 

other fields is very simple: at n = there arc no mixings in (6.6) and the only source of corrections is that T^^ is 
coupled to shifted h^i, + jz2^??;lii/ instead of T^i^. This provides an additional contribution from the 0-exchange, 
which is equal to 



I V d-2T^ 1 



d-2j d-lfc2 {d-l){d-2)k^ 



(6.14) 



(the second factor takes into account the coefficient in front of kinetic term for cj)) and should be added to the 
pure graviton exchange (6.13), thus changing —335 ^° + {d-i)(d-2) ~ which reproduces (6.12). 

In other words, we observe that in the KK theory the current- current interaction gets contributions from 3 
channels, 

KK theory = graviton + ghost + scalar field (6.15) 
while in the massless gravity these are 2 channels 



Massless gravity = graviton + ghost (6.16) 

and in the PF theory, where the ghost is removed out of the spectrum by taking its mass to infinity, there is 
only the graviton contribution: 

PF theory = graviton (6-17) 

As we saw above, the contributions of the ghost and the scalar field exactly cancel each other so that the 
current- current interactions in the KK theory and in the PF theory coincide. If, however, one considers a KK 
theory with several compactificd dimensions, i.e. with several scalar fields added, this compensation will no 
longer take place, and all the three interactions, (6.15), (6.16) and (6.17) will be different. 
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7 Massive gravity within and beyond the quadratic approximation 



In this section we make very brief comments about the other aspects of massive gravity, some of which are not 

directly seen at linearized level, but are in fact direct consequences of the properties of quadratic theory. These 
are pronounced and affective phenomena, but do not add much new to the theoretical aspects of the problem. 

7.1 Vainshtein radius [8] 

This celebrated result is often considered as the clear proof of pathology of PF massive gravity, though the actual 
statement is much simpler. The new modes, revived by addition of mass terms to Einstein-Hilbert action, get 
kinetic terms because at least some of them are Stueckelberg fields. We emphasize once again that one does 
not modify the kinetic part of Einstein-Hilbert action, only masses (more generally, a non-trivial potential) are 
added. This means that the new kinetic terms enter with small coefficients m^, m^fc^s^, and this means that 
the physical fields are ms rather than s. This means that higher-degree terms, say, m^s"^ = ^{ms)^ are actually 
entering with large couplings m~^, and that the theory is actually strongly coupled, moreover interactions 
get stronger with decrease of m^. Interaction effects arc of course falling with the distance, therefore very far 
away, beyond some "Vainshtein radius" perturbative regime can still be reliable, but in general perturbation 
theory does not make much sense at small and even moderate distances. 

This strong coupling effect is a result of explicit breakdown of gauge symmetry and it is rather similar 
to what happens when Higgs mass goes beyond the unitary limit in the case of spontaneously broken gauge 
symmetry in Yang-Mills theory. Nothing like this strong-coupling phenomena occurs in the case of gauge 
invariant Kaluza-Klein massive gravity and hence it can be avoided when gravity mixes with other fields. 

7.2 Ghosts 

Ghost are the fields which enter Lagrangian with the wrong sign in front of the time-derivative term (jp. What 
happens to them is that they grow in time until this growth is stopped by non-linear terms in the action or the 
theory gets into a different phase with another spectrum of quasiparticles. In this sense this is a phenomenon of 
the same class as the previous one: ghosts make perturbative treatment unreliable and most often inadequate. 
In this context one often speaks about "negative norms" and "violation of unitarity" - but this actually refers to 
" perturbative unitarity" , implying that the actual spectrum of the theory is more-or-less accurately described 
by quadratic part of the Lagrangian, what is not always the case. While occurrence of ghosts clearly means 
that the theory is not what we thought it will be, it does not obligatory means that it is ill and un-curable - it 
is just not perturbative and most probably strongly coupled. 

There are three ways to deal with ghosts. The first option is to study the theory as it is (what is rarely 
done, see, however. [31] or the series of works about the "pathological" Hamiltonian H = xp [32]). The second 
option is to say that the growth of ghosts is slow enough to be acceptable (e.g., the ghost formally appears 
in the Lorentz invariant massive gravity with A = 2B, which corresponds to addition of the cosmological 
term to the Einstein-Hilbert Lagrangian, what happens is that the fields grow together with the growth of the 
Universe itself). The third, and most popular option is to fine-tune parameters in order to "eliminate" ghosts, 
for example, to give them an infinite mass. 

7.3 BoulwEire-Deser mode [9, 13] 

The last option is exactly the one, "distinguishing" the Pauli-Fierz gravity with A = B. However, the way in 
which the ghost acquires infinite mass is somewhat special and actually unreliable. The mass is infinite because 
the coefficient in front of kinetic term vanishes! As we already discussed, this is what actually makes the theory 
strongly coupled, moreover, this of course makes the fine-tuning fully unreliable. Any minor deformation of the 
theory destroys the fine-tuning and brings the ghost back to existence. Boulware-Deser instability is a concrete 
example of this phenomenon: switching on non-trivial background metric contributes to the coeflacient in front 
of the kinetic term and shifts it away from zero. Since background is arbitrary it can not be compensated by 
variation of just two adjusting constants A and B. 

Of course, nothing like this happens in ghost-free generalizations of PF gravity such as KK or other models 
with extra fields. Ghosts are absorbed into these additional fields in a universal, background independent way. 

7.4 DVZ discontinuities 

As we already explained these discontinuities occur when comparing two theories with different sets of fields and 
thus are of no surprise. However, from the point of view of concrete example of massless gravity, it is instructive 
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to distinguish between two situations: when one compares it to theory with more fields and with less fields. 

The first is, for example, the ease of Kaluza-Klein gravity: in a given mass sector of Kaluza-Klcin theory 
there are two more fields, and 0, which mix with the massive graviton and also contribute to Newton-like 
interaction. This additional contribution explains the difiference. 

The second is the Pauli-Fierz gravity, where one of the modes, contributing to Newton interaction in massless 
theory, becomes a ghost and is thrown away " by hands" . This explains the discontinuity. It is in no way a 
property of massive gravity, it is a property of this artificial throw-away prescription. The reason why the 
would-be ghost is allowed in massless gravity is that it has the same mass as the other (non-ghost) degrees 
of freedom, and this results into decoupling of the negative norm states from the spectrum which leaves the 
theory perturbatively consistent (we do not speak about UV problems of quantum gravity here). Notably, the 
same kind of absorption could be expected if instead of the Pauli-Fierz choice A = B we put A = 2B when 
the two masses continue to coincide, M = m. Then there will be no DVZ discontinuity. This is the case of 
the cosmological term added to the Einstein-Hilbert action, and, because of presence of linear terms in the 
Lagrangian, one has to re-expand this latter above the non-flat AdS vacuum where the linear term cancels [33] 
(see also footnote 7 in [13] about this option). 

7.5 Tachyons 

The name tachyon refers to superluminal propagation. However, in modern literature it is actually used to 
mean something different: a mode, that reflects perturbative instability of the backgromid around which the 
perturbative theory is developed. Technically this means that there is a pole in the propagator at vanishing 
frequency, for example, < 0. Physically this means that a phase transition of the first kind occurs from 
perturbatively unstable vacuum to another one, stable at least perturbatively (false vacua which are separated 
from the real ones by potential barriers do not have tachyons in perturbative spectra, their instability is essen- 
tially non-perturbative phenomenon). Such phase transition takes place spontaneously and independently in 
all point of the space: this can look like a propagation of non-causal (superluminal) particle, but has a clear 
reason, and processes taking place in different places are indeed casually unrelated. Tachyons are in no way a 
problem of the theory they signal just that we treat (interpret) it in a wrong way. It is of course not always 
easy to find the write vacuum, this is often related to finding the write non-perturbative formulation of the 
theory. The celebrated example of such lasting study was interpretation of tachyons in string theory (for open 
superstring the puzzle was partly resolved by A. Sen in [34], in general string models it remains obscure). 

What is the right vacuum of Lorentz- violating massive gravity with tachyons, and what at all is its adequate 
non-perturbative formulation is an open problem (not surprisingly given by the young age of this subject). 
However, one may expect this vacuum has not to be homogeneous, due to non-trivial dispersion laws (see s.5.6). 
Moreover, phenomenologically this may be not that bad, since the scale of these inhomogeneities is determined 
by inverse masses that breaks Lorentz invariance, i.e. which are very small. Inhomogeneities at such large scales 
are quite possible. 

7.6 Superluminal propagation [26, 23] 

We reserve the word superluminal propagation for phenomenon which is (at least looks) different from tachyons, 
i.e. is not related in any clear way to perturbative instabilities and phase transitions. This is occurrence of 
poles in the propagator at uj = ck with c > 1. Such poles are forbidden by Lorentz-invariance, but we saw 
that they naturally appear in the spectrum of massive gravity when Lorentz invariance is violated (for example, 
= mf/m2). The meaning of superluminal dispersion laws remains controversial, we feel that most people 
find them unpleasant, referring mostly to causality arguments, which are also used against time-machines - 
which surely are not forbidden by the laws of nature (see [22] for the most recent discussion). Quite similarly, 
super luminals are unavoidably present in our theories, whether we like them or not. Whenever one perturbs a 
light-like dispersion law ut^k'^ for a collection of light-like particles (two polarizations are already sufficient) in 
a Lorentz- violating way, the eigenvalues split, and one goes above, another below c = 1: this is a fundamental 
law of linear algebra, well known in the level-splitting theory in quantum mechanics. It is important that 
Lorentz violation need not be "fundamental" (i.e. explicitly written down into Lagrangian): a Lorentz- violating 
background is already enough. The celebrated example is emergency of superluminal photons in curved space 
[26], where the effective Lagrangian acquires a quantum correction (1-loop of fermion or any other field in 
gravitational background) and becomes 




(7.1) 
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what - according to the above-mentioned Unear-algebra theorem - unavoidably leads to dispersion law with 

c > 1 for one of photon polarizations. 

Implications of supcrluminals arc still badly understood. It is not oven fully clear if they can be used to 
construct time-machines, either big (of astrophysical scale) [35] or mini (of Planckian scale) [22]. In any case, 
we repeat that time machines are allowed already in ordinary General Relativity and, within hypothetical TeV- 
gravity models [20], might even be massively created (though immediately evaporate) in particle collisions in 
modern accelerators [22]. 

Most important, our intuition about Lorentz-non-invariant theories is still very underdeveloped. Always 
we discuss our models (in this context) from the point of view of some outside observer which believes into 
fundamental Lorentz invariance: we ask what happens if we look on a superluminal from another frame (what 
happens is a singularity in dispersion rule uj + [3k = c{k + /Jlo) i.e. uj = at /3 = 1/c for c > 1 instead of 

the usual zero at /? = c for c < 1) and discuss whether this good or bad. However, if we just look at a theory 
per se, it is absolutely unclear what is going to be bad about it: there simply is no Lorentz invariance and it is 
unclear why at all one should ask what happens if you make such a transformation. It is like writing down a 
theory with rescaled a;-coordinate and apply Lorentz transformation with and old (not-rescaled) x. Even more 
important would be to examine what can be wrong with just a combination of two fields, 

L = li^l - clV4>l +4>l- C2W2) + V{<ki,4>2) (7.2) 

If C2 > ci the second field is a superluminal from the point of view of the first one, but nothing wrong is expected 
if we look on the first field from the point of view of the second. 

7.7 Radiation of massive gravitons 

As mentioned in the introduction, one of our original questions was if the "pathologies" of massive gravity 

can somehow affect the analysis [19] of massivc-graviton radiation in TeV-gravity models, a worry, naturally 
implied by the original analysis in [18]. However, the Kaluza-Klein theory (which is in the base of the Tev- 
gravity models) seems to be free of any problems, true or imaginary, of generic massive gravity and our current 
feeling is that one can treat radiation of massive Kaluza-Klein gravitons in the straightforward and naive way, 
as suggested [19]. The problem, however, deserves an independent analysis. 
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